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Abstract 

We provide a systematic and self-contained exposition of the subject of localized qubits in curved 
spacetimes. This research was motivated by a simple experimental question: if we move a spatially 
localized qubit, initially in a state along some spacetime path T from a spacetime point x\ to 

another point xi, what will the final quantum state \if>-i) be at point £2? This paper addresses this 
question for two physical realizations of the qubit: spin of a massive fermion and polarization of a 
photon. Our starting point is the Dirac and Maxwell equations that describe respectively the one- 
particle states of localized massive fermions and photons. In the WKB limit we show how one can 
isolate a two-dimensional quantum state which evolves unitarily along Y. The quantum states for these 
two realizations are represented by a left-handed 2-spinor in the case of massive fermions and a four- 
component complex polarization vector in the case of photons. In addition we show how to obtain from 
this WKB approach a fully general relativistic description of gravitationally induced phases. We use this 
formalism to describe the gravitational shift in the Colella-Overhauser-Werner 1975 experiment. In the 
non-relativistic weak field limit our result reduces to the standard formula in the original paper. We 
provide a concrete physical model for a Stern-Gerlach measurement of spin and obtain a unique spin 
operator which can be determined given the orientation and velocity of the Stern-Gerlach device and 
velocity of the massive fermion. Finally, we consider multipartite states and generalize the formalism to 
incorporate basic elements from quantum information theory such as quantum entanglement, quantum 
teleportation, and identical particles. The resulting formalism provides a basis for exploring precision 
quantum measurements of the gravitational field using techniques from quantum information theory. 
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Notation and conventions 

We use the following index notation: 

- fj,, v, p, a, . . . denote spacetime tensor indices 

- I, J,K, L, . . . = 0, 1, 2, 3 denote tetrad indices. 

- i,j, k,l, . . . = 1, 2, 3 for spatial components of the tetrad (the 'triad') 

- A, B, C, D, . . , = 1, 2 for spinor indices 

- A', B', C, D', . . . = 1, 2 for conjugate spinor indices 

The Minkowski metric is defined as rj^ — diag(l, — 1, —1, —1). We generally use natural units where c 
h = 1, and in addition we set the charge of a proton to e = 1. 
We use the Weyl representation for the Dirac 7-matrices 

1 ( a 1 

where a 1 = and a 1 = (1, — u l ), and a 1 are the usual Pauli matrices. Writing this object in spinor 

notation we have a 1 = a 1 AA , and a 1 = a IA . In order to interpret the spatial parts a 1 AA , and a %A A as 
the Pauli matrices we use the convention in pQ: for a 1 = a IAA the primed index is the row index and 
the unprimed index is the column index, and the opposite assignment occurs for a 1 = AA ,. In spinorial 
notation a 1 is not an operator, rather an operator A carries an index structure A A or A A B . Throughout 
this paper we will switch between the implicit index notation A and A A or A A B . 



1 Introduction 

This paper will provide a systematic and self-contained exposition of the subject of localized qubits in 
curved spacetimes with the focus on two physical realizations of the qubit: spin of a massive fermion and 
polarization of a photon. Although a great amount of research has been devoted to quantum field theory in 
curved spacetimes G3 0] and also more recently to relativistic quantum information theory in the presence 
of particle creation and the Unruh effect El El HI [9] , the literature about localized qubits and quantum 
information theory in curved spacetimes is relatively sparse |10 [ 111 ) H"2"]. In particular, we are aware of only 
three papers, [101 EH US], that deal with the following question: if we move a spatially localized qubit, 
initially in a state l^i), along some spacetime path T from a point p\ in spacetime to another point P2, what 
will the final quantum state \tjj2} be at point P2? This, and other relevant questions, were given as open 
problems in the field of relativistic quantum information by Peres and Terno in [141 p. 19]. The formalism 
developed in this paper will be able to address such questions, and will also be able to deal with the basic 
elements of quantum information theory such as entanglement and multipartite states, teleportation, and 
quantum interference. 

The basic object in quantum information theory is the qubit. Given a Hilbert space of some physical 
system, we can physically realize a qubit as any two-dimensional subspace of that Hilbert space. However, 
such physical realizations will in general not be localized in physical space. We shall restrict our attention 
to physical realizations that are well-localized in physical space so that we can approximately represent the 
qubit as a two-dimensional quantum state attached to a single point in space. From a spacetime perspective 
a localized qubit is then mathematically represented as a sequence of two-dimensional quantum states along 
some spacetime trajectory corresponding to the worldline of the qubit. 

In order to ensure relativistic invariance it is then necessary to understand how this quantum state 
transforms under a Lorcntz transformation. However, as is well-known, there are no finite-dimensional 
faithful unitary representations of the Lorentz group [T5] and in particular no two-dimensional ones. The 
only faithful unitary representations of the Lorentz group are infinite dimensional (see e.g. [16 ). Hence, 
these cannot be taken to mathematically represent a qubit, i.e. a two-level system. Naively it would appear 
that a formalism for describing localized qubits which is both relativistic and unitary is a mathematical 
impossibility. 



4 



In the case of flat spacetime the Wigner representations [TH [T7] provide unitary and faithful but infinite- 
dimensional representations of the Lorentz group. These representations make use of the symmetries of 
Minkowski spacetime, i.e. the full inhomogeneous Poincare group which includes rotations, boosts, and 
translations. The basis states \p, a) are taken to be eigenstates of the four momentum operators (the gener- 
ators of spatio-temporal translations) P M , i.e. P^\p, a) = p^lp, a) where the symbol a refers to some discrete 
degree of freedom, perhaps spin or polarization. One strategy for obtaining a two-dimensional (perhaps 
mixed) quantum state p ac! i for the discrete degree of freedom a would be to trace out the momentum degree 
of freedom. But as shown in [T5J []1|1 EH HO] this density operator does not have covariant transformation 
properties. The mathematical reason, from the theory presented in this paper, is that the quantum states 
for qubits with different momenta belong to different Hilbert spaces. Thus, the density operator p aa i is then 
a mixture of states which belong to different Hilbert spaces. The operation of 'tracing out the momenta' is 
neither physically meaningful nor mathematically motivated. 

Another strategy for defining qubits in a relativistic setting would be to restrict to momentum eigenstates 
\p, a). The continuous degree of freedom P is then fixed and the remaining degrees of freedom arc discrete. 
In the case of a photon or fermion the state space is two dimensional and this can then serve as a relativistic 
realization of a qubit. This is the strategy in [TUHH] where the authors develop a theory of transport of qubits 
along worldlines. However, when we go from a flat spacetime to curved we lose the translational symmetry 
and thereby also the momentum eigenstates \p, a) . The only symmetry remaining is local Lorentz invariance 
which is manifest in the tetrad formulation of general relativity. Since the translational symmetry is absent 
in a curved spacetime it seems difficult to work with Wigner representations which rely heavily on the full 
inhomogeneous Poincare group. The use of Wigner representations therefore needs further justification as 
they do not exist in curved spacetimcs. 

In this paper we shall refrain altogether from making use of the infinite-dimensional Wigner representa- 
tion. Since our focus is on qubits physically realized as polarization of photons and spin of massive fermions 
our starting point will be the field equations that describe those physical systems, i.e. the Maxwell and 
Dirac equations in curved spacetimes. Using the WKB approximation we then show in detail how one can 
isolate a two-dimensional Hilbert space and determine an inner product, unitary evolution, and a quantum 
state. Our procedure reproduces the results of [101 111] , and can be regarded as an independent justification 
and validation. 

Notably, possible gravitationally induced global phases [21, 22, 23 , 24, 25, 26 , which are absent in [TUllllj . 
are automatically included in the WKB approach. Such a phase is irrelevant if only single trajectories are 
considered. However, quantum mechanics allows for more exotic scenarios such as when a single qubit is 
simultaneously transported along a superposition of paths. In order to analyze such scenarios it is necessary 
to determine the gravitationally induced phase difference. We show how to derive a simple but fully general 
relativistic expression for such a phase difference in the case of spacetime Mach-Zehnder interferometry. Such 
a phase difference can be measured empirically [27] with neutrons in a gravitational field. See (23] HH1 US] 
and references therein for further details and generalizations. The formalism developed in this paper can 
easily be applied to any spacetime, e.g. spacetimes with frame-dragging. 

This paper aims to be self-contained and we have therefore included necessary background material such 
as the tetrad formulation of general relativity, the connection 1-form, spinor formalism and more (see |I]and 
[A| . For example, the absence of global reference frames in a curved spacetime has a direct bearing on how 
entangled states and quantum teleportation in a curved spacetime are to be understood conceptually and 
mathematically. We discuss this in section [9] 

2 An outline of methods and concepts 

In this section we provide a general outline of the main ideas and concepts needed to understand the topic 
of localized qubits in curved spacetimes. 

2.1 Localized qubits in curved spacetimes 

Let us now make precise the concept of a localized qubit. As a minimal characterization, a localized qubit 
is understood in this paper as any two-level quantum system which is spatially well-localized. Such a qubit 
is effectively described by a two-dimensional quantum state attached to a single point in space. From a 
spacetime perspective the history of the localized qubit is then a sequence (i.e. a one-parameter family) 
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of two-dimensional quantum states each associated with a point a; M (A) on the worldline of the qubit 

parameterized by A. In this paper we will focus on qubits represented by the spin of an electron and the 
polarization of a photon and show how one can, by applying the WKB approximation to the corresponding 
field equation (the Dirac or Maxwell equation) , extract a two- level quantum state associated with a spatially 
localized particle. 

The sequence of quantum states \tp(X)) must be thought of as belonging to distinct Hilbert spaces 1-L X (X) 
attached to each point a; M (A) of our trajectory. The situation is identical to that in differential geometry where 
one must think of the tangent spaces associated with different spacetime points as mathematically distinct: 
since the parallel transport of a vector along some path from one point to another is path dependent there is 
no natural identification between vectors of one tangent space and the other. The parallel transport, for any 
type of object, is simply a sequence of infinitesimal Lorentz transformations acting on the object and it is 
this sequence that is in general path dependent. Thus, if we are dealing with a physical realization of a qubit 
whose state transforms non-trivially under the Lorentz group, as is the case for the two physical realizations 
that we are considering, we must also conclude that in general it is not possible to compare quantum states 
associated with distinct points in spacetime. As we shall see in sections 5.3.1 and 6.5 Hilbert spaces for 
different momenta p^{\) of the particle carrying the qubit must also be considered distinct. The Hilbert 
spaces will therefore be indexed as H X n, and so along a trajectory there will be a family of Hilbert spaces 

The ambiguity in comparing separated states has particular consequences: It is in general not well- 
defined to say that two quantum states associated with distinct points in spacetime are the same. Nor is 
it mathematically well-defined to ask how much a quantum state has "really" changed when moved along a 
path. Nevertheless, if two initially identical states are transported to some point x but along two distinct 
paths, the difference between the two resulting states is well-defined, since we are comparing states belonging 



to the same Hilbert space (see figure 2.1 1. 

There are also consequences for how we interpret basic quantum information tasks such as quantum 
teleportation: When Alice "teleports" a quantum state over some distance to Bob we would like to say that 
it is the same state that appears at Bob's location. However, this will not have an unambiguous meaning. An 
interesting alternative is to instead use the maximally entangled state to define what is "the same" quantum 
state for Bob and Alice, at their distinct locations. We return to these issues in j]9.3| 




Figure 2.1: If we parallel transport a vector v from point x± to X2 along two distinct trajectories 
Ti and F2 in curved spacetime we generally obtain two distinct vectors vr 1 and vr 2 at X2- Thus, 
no natural identification of vectors of one tangent space and another exists in general and we need 
to associate a distinct tangent space for each point in spacetime. The same applies to quantum 
states and their Hilbert spaces: the state at a point X2 of a qubit moved from a point x\ would 
in general depend on the path taken, and hence the Hilbert space for each point is distinct. The 
Hilbert spaces H X1 and H X2 are illustrated as vertical 'fibres' attached to the spacetime points 
xi and X2- 

In a strict sense a localized qubit can be understood as a sequence of quantum states attached to 
points along a worldline. We will however relax this notion of localized qubits slightly to allow for path 
superpositions as well. More specifically, we can consider scenarios in which a single localized qubit is split 
up into a spatial superposition, transported simultaneously along two or more distinct worldlines, and made 
to recombine at some future spacetime region so as to produce quantum interference phenomena. We will 
still regard these spatial superpositions as localized if the components of the superposition are each localized 
around well-defined spacetime trajectories. 
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2.2 Physical realizations of localized qubits 

The concepts of a classical bit and a quantum bit (cbit and qubit for short) are abstract concepts in the 
sense that no importance is usually attached to the specific way in which we physically realize the cbit or 
qubit. However, when we want to manipulate the state of the cbit or qubit using external fields, the specific 
physical realization of the bit becomes important. For example, the state of a qubit, physically realized as 
the spin of a massive fermion, can readily be manipulated using an external electromagnetic field, but the 
same is not true for a qubit physically realized as the polarization of a photon. 

The situation is no different when the external field is the gravitational one. In order to develop a 
formalism for describing transport of qubits in curved spacetimes it is necessary to pay attention to how the 
qubit is physically realized. Without knowing whether the qubit is physically realized as the spin of a massive 
fermion or the polarization of a photon, for example, it is not possible to determine how the quantum state 
of the qubit responds to the gravitational field. More precisely: gravity, in part, acts on a localized qubit 
through a sequence of Lorentz transformations which can be determined from the trajectory along which 
it is transported and the gravitational field, i.e. the connection one-form oj^ 1 j. Since different qubits can 
constitute different representations under the Lorentz group, the influence of gravity will be representation 
dependent. This is not at odds with the equivalence principle, which only requires that the qubits are acted 
upon with the same Lorentz transformation. 



2.3 Our approach 

Our starting point will be the one-particle excitations of the respective quantum fields. These one-particle 
excitations are fields "J or A^, which are governed by the classical Dirac or Maxwell equation, respectively. 
Our goal is to formulate a mathematical description for localized qubits in curved spacetime. Therefore we 
must find a regime in which the spatial degrees of freedom of the fields are suppressed so that the relevant 
state space reduces to a two-dimensional quantum state associated with points along some well-defined 
spacetime trajectory. Our approach is to apply the WKB approximation to these field equations (sections 



5.1 and 6.1) and study spatially localized solutions. In this way we can isolate a two-dimensional quantum 
state that travels along a classical trajectory. 

In the approach that we use for the two realizations, we start with a general wavefunction for the fields 
expressed as 

cj> A (x) = ^ A (x)v(x)e i9 ^ or A^x) = Re[^(xMx)e ie ^} 

where the two-component spinor field 4>a(x) is the left-handed component of the Dirac field \Jj and x^ 
is some coordinate system. The decompositions for the two fields are similar: 9{x) is the phase, ip(x) is 
the real- valued envelope, and iPa(x) or ipfj.i x ) are fields that encode the quantum state of the qubit in 
the respective cases. These latter objects are respectively the normalized two-component spinor field and 
normalized complex-valued polarization vector field. Note that we are deliberately using the same symbol ip 
for both the two-component spinor tp A and the polarization 4-vector ip^ as it is these variables that encode 
the quantum state in each case. 

The WKB limit proceeds under the assumptions that the phase 0(x) varies in x much more rapidly 
than any other aspect of the field and that the wavelength of the phase oscillation is much smaller than the 
spacetime curvature scale. Expanding the field equations under these conditions we obtain: 

• a field of wavevectors k^(x) whose integral curves satisfy the corresponding classical equations of 
motion; 

• a global phase 6, determined by integrating fc p along the integral curves; 

• transport equations that govern the evolution of ipA and ip^ along this family of integral curves; 

• a conserved current which will be interpreted as a quantum probability current. 

The assumptions of the WKB limit by themselves do not ensure a spatially localized envelope (p(x), and 



therefore do not in general describe localized qubits. In sections 5.2 and 6.2 we add further assumptions 
that guarantee that the qubit is localized during its transport along the trajectory. The spatial degrees 
of freedom are in this way suppressed and we can effectively describe the qubit as a sequence of quantum 
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states, encoded in the objects ipA ( T ) or VvM- These objects constitute non-unitary representations of the 
Lorentz group. As we shall see in Sj5J unitarity is recovered once we have correctly identified the respective 
inner products. Notably, the Hilbert spaces T-L( x ,p)(X) we obtain are labelled with both the position and the 
momentum of the localised qubit. 

Finally, since the objects ?Pa(t) and VvM hare been separated from the phase e lB ^ x \ the transport 
equations for these objects do not account for possible gravitationally induced global phases. We show how 
to obtain such phases in Sj7] from the WKB approximation. Thus, with the inclusion of phases, we have 
provided a complete, Lorentz covariant formalism describing the transport of qubits in curved spacetimes. 
Hereafter it is straightforward to extend the formalism to several qubits in order to treat multipartite states, 
entanglement and teleportation (Q, providing the basic ingredients of quantum information theory in curved 
spacetimes. 



3 Issues from quantum field theory and the domain of applicabil- 
ity 

The formalism describing qubits in curved spacetimes presented in this paper has its specific domain of 
applicability and cannot be taken to be empirically correct in all situations. One simple reason for this is 
that the current most fundamental theory of nature is not formulated in terms of localized qubits but instead 
involves very different objects such as quantum fields. There are four important issues arising from quantum 
field theory that restrict the domain of applicability: 

• the problem of localization; 

• particle number ambiguity; 

• particle creation; 

• the Unruh effect. 

Below we discuss these issues and indicate how they restrict the domain of applicability of the formalism of 
this paper. 



3.1 The localization problem 

The formalism of this paper concerns spatially localized qubits, with the wavepacket width being much 
smaller than the curvature scale. However, it is well-known from quantum field theory that it is not possible 
to localize one-particle states to an arbitrary degree. For example, localization of massive fermions is limited 
by the Compton wavelength A c = h/mc 30 . More precisely, any wavefunction constructed from exclusively 
positive frequency modes must have a tail that falls off with radius r slower than e _r / Ac . However, this 
is of no concern if we only consider wavepackets with a width much larger than the Compton wavelength. 
This consequently restricts the domain of applicability of the material in this paper. In particular, since the 



width of the wavepacket is assumed to be much smaller than the curvature scale (see [3.2), the localization 
theorem means that we cannot deal with extreme curvature scales of the order of the Compton wavelength. 

A similar problem exists also for photons. Although the Compton wavelength for photons is ill-defined, 
it has also been shown that they must have non- vanishing sub-exponential tails [31] [32] • 

Given these localization theorems it is not strictly speaking possible to define a localized wavepacket with 
compact support. However, for the purpose of this paper we will assume that most of the wavepacket is 
contained within some region, smaller than the curvature scale, and the exponential tails outside can safely 
be neglected in calculations. We will assume from here on that this is indeed the case. 



3.2 Particle number ambiguity 

One important lesson that we have learned from quantum field theory in curved spacetimes is that a natural 
notion of particle number is in general absent; see e.g. [3]. It is only under special conditions that a natural 
notion of particle number emerges. Therefore, for arbitrary time-dependent spacetimes it is not in general 
possible to talk unambiguously about the spin of one electron or the polarization state of one photon as this 
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would require an unambiguous notion of particle number. This is important in this paper because a qubit 
is realized by the spin of one massive fermion or polarization of one photon. 

The particle number ambiguity can be traced back to the fact that the most fundamental mathematical 
objects in quantum field theory are the quantum field operators and not particles or Fock space representa- 
tions. More specifically, how many particles a certain quantum state is taken to represent depends in general 
on how we expand the quantum field operators in terms of annihilation and creation operators (d^, aj): 

4>( x ) = ^2 fi&i + fi^i 
i 

which in turn depends on how the complete set of modes (which are solutions to the corresponding classical 
field equations) is partitioned into positive and negative frequency modes (ft, fi). In particular, the number 
operator N = J^i a l^i depends on the expansion of the quantum field operator (f>(x). This expansion can be 
done in an infinitude of distinct ways related by Bogoliubov transformations [2] . Particle number is therefore 
ill-defined. Since we base our approach on the existence of well-defined one-particle states for photons and 
massive fcrmions, the particle number ambiguity seems to raise conceptual difficulties. 

We will now argue from the equivalence principle that the particle number ambiguity does not occur for 
spatially localized states. Consider first vanishing external fields and thus geodesic motion (we will turn to 
non-geodesics in the next section). In a pseudo-Riemannian geometry, for any sufficiently small spacetime 
region we can always find coordinates such that the metric tensor is the Minkowski metric g^ v = r\av and 
the affine connection is zero TZ V = 0. However, this is true also for a sufficiently narrow strip around any 
extended spacetime trajectory, i.e. there exists an extended open region containing the trajectory such that 
9iiv — f)nv and = [33] . Thus, as long as the qubit wavepacket is confined to that strip it might as 
well be travelling in a flat spacetime. In fact, the usual free Minkowski modes e ±zp ' x form a complete set of 
solutions to the wave equation for wavepackets localized within that strip. Using these modes we can then 
define positive and negative frequency and thus the notion of particle number becomes well-defined. Thus, if 
we restrict ourselves to qubit wavepackets that are small with respect to the typical length scale associated 
with the spacetime curvature, the particle number ambiguity is circumvented and it becomes unproblematic 
to think of the classical fields ^(x) and A^{x) as describing one-particle excitations of the corresponding 
quantum field. 

3.3 Particle creation and external fields 

Within a strip as defined in the previous section, the effects of gravity are absent and therefore there is no 
particle creation due to gravitational effects for sufficiently localized qubits. If the trajectory T along which 
the qubit is transported is non-geodesic, non-zero external fields need to be present along the trajectory. 
For charged fermions we could use an electromagnetic field. However, if the field strength is strong enough 
it might cause spontaneous particle creation and we would not be dealing with a single particle and thus 
not a two-dimensional Hilbert space. As the formalism of this paper presupposes a two-dimensional Hilbert 
space, we need to make sure that we are outside the regime where particle creation can occur. 

When time-dependent external fields are present, the normal modes e ±ip ' x are no longer solutions of the 
corresponding classical field equations and there will in general be no preferred way of partitioning the modes 
(/i: fi) m t° positive and negative frequency modes. Therefore, even when we confine ourselves to within the 
above mentioned narrow strip, particle number is ambiguous. 

This type of particle number ambiguity can be circumvented with the help of asymptotic 'in' and 'out' 
regions in which the external field is assumed to be weak. In the scenarios considered in this paper there will 
be a spacetime region TZ prep . in which the quantum state of the qubit is prepared, and a spacetime region 
TZ-meas. where a suitable measurement is carried out on the qubit. The regions are connected by one or many 
timelike paths along which the qubit is transported. The are here taken to be 

macroscopic but still sufficiently small such that no tidal effects are detectable, and so special relativity is 
applicable. We allow for non-zero external fields in these regions and along the trajectory, though we assume 
that external fields (or other interactions) are weak in these end regions so that the qubit is essentially free 
there. This means that in TZ prep . and lZ meas . we can use the ordinary Minkowski modes e lpx and e~ lp ' x 
to expand our quantum field. This provides us with a natural partitioning of the modes into positive and 
negative frequency modes and thus particle number is well-defined in the two regions TZprep. an d 7\L meas .. For 
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our purposes we can therefore regard (approximately) the regions TZ pre p. and TZ-meas. as the asymptotic 'in' 
and 'out' regions of ordinary quantum field theory. 

If we want to determine whether there is particle creation we simply 'propagate' (using the wave equation 
with an external field) a positive frequency mode (with respect to the free Minkowski modes in lZ pre p.) from 
region lZ pre p. to 1Z me as.- In region lZ me as. we then see whether the propagated mode has any negative fre- 
quency components (with respect to the free Minkowski modes in lZ meas ,). If negative frequency components 
are present we can conclude that particle creation has occurred (see e.g. [34] ) . This will push the physics 
outside our one-particle-excitation formalism and we need to make sure that the strength of the external 
field is sufficiently small so as to avoid particle creation. 

One also has to avoid spin-flip transitions in photon radiation processes such as gyromagnetic emission, 
which describes radiation due to the acceleration of a charged particle by an external magnetic field, and 
the related Brcmsstrahlung, which corresponds to radiation due to scattering off an external electric field 
[331 1551 157] . For the former, a charged fermion will emit photons for sufficiently large accelerations and can 
cause a spin flip and thus a change of the quantum state of the qubit. Fortunately, the probability of a 
spin-flip transition is much smaller than that of a spin conserving one, which docs not alter the quantum 
state of the qubit [37] . In this paper we assume that the acceleration of the qubit is sufficiently small so that 
we can ignore such spin-flip processes. 

3.4 The Unruh effect 

Consider the case of flat spacetime. A violently accelerated particle detector could click (i.e. indicate that 
it has detected a particle) even though the quantum field <p is in its vacuum state. This is the well-known 
Unruh effect [35] [U [5] . What happens from a quantum field theory point of view is that the term for the 
interaction between a detector and a quantum field allows for a process where the detector gets excited 
and simultaneously excites the quantum field. This effect is similar to that when an accelerated electron 
excites the electromagnetic field [39]. A different way of understanding the Unruh effect is by recognizing 
that there are two different timelike Killing vector fields of the Minkowski spacetime: one generates inertial 
timelike trajectories and the other generates orbits of constant proper acceleration. Through the separation of 
variables of the wave equation one then obtains two distinct complete sets of orthonormal modes: Minkowski 
modes and Rindler modes, corresponding respectively to each Killing field. The positive Minkowski modes 
have negative frequency components with respect to the Rindler modes and it can be shown that the 
Minkowski vacuum contains a thermal spectrum with respect to a Rindler observer. 

In order to ensure that our measurement and preparation devices operate 'accurately', their acceleration 
must be small enough so as not to cause an Unruh type effect. 

3.5 The domain of applicability 

Let us summarize. In order to avoid unwanted effects from quantum field theory we have to restrict ourselves 
to scenarios in which: 

• the qubit wavepacket size is much smaller than the typical curvature scale (to ensure no particle number 
ambiguity) ; 

• in the case of massive fermions, because of the localization problem the curvature scale must be much 
larger than the Compton wavelength; 

• there is at most moderate proper acceleration of the qubit (to ensure no particle creation or spin-flip 
transition due to external fields); 

• there is at most moderate acceleration of preparation and measurement devices (to ensure negligible 
Unruh effect). 

For the rest of the paper we will tacitly assume that these conditions are met. 

4 Reference frames and connection 1-forms 

The notion of a local reference frame, which is mathematically represented by a tetrad field ej(x), is essential 
for describing localized qubits in curved spacetimes. This section provides an introduction to the mathematics 
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of tetrads with an eye towards its use for quantum information theory in curved spacetime. The hurried 
reader may want to skip to fj5] A presentation of tetrads can also be found in jJDl App. J]. 

4.1 The absence of global reference frames 

One main issue that arises when generalizing quantum information theory from flat to curved spaces is the 
absence of a global reference frame. On a flat space manifold one can define a global reference frame by first 
introducing, at an arbitrary point X\, some orthonormal reference frame, i.e. we associate three orthonormal 
spatial vectors (x Xl , y Xl , z Xl ) with the point x\. In order to establish a reference frame at some other point 
X2 we can parallel transport each of the three vectors to that point. Since the manifold is flat the three 
resulting orthonormal directions are independent of the path along which they were transported. Repeating 
this for all points x in our space we obtain a unique field of reference frames (x x , y x ,z x ) defined for all points 
x on the manifoldf^] Thus, from an arbitrarily chosen reference frame at a single point X\ we can erect a 
unique global reference frame. 

However, when the manifold is curved no unique global reference frame can be established in this way. 
The reference frame obtained at point xi by the parallel transport of the reference frame at x\ is in general 
dependent on the path along which the frame was transported. Thus, in general there is no path-independent 
way of constructing global reference frames. Instead we have to accept that the choice of reference frame at 
each point on the manifold is completely arbitrary, leading us to the notion of local reference frames. 

To illustrate this situation and its consequences in the context of quantum information theory in curved 
space, consider two parties, Alice and Bob, at separated locations. First we turn to the case where the 
space is flat and the entangled state is the singlet state. The measurement outcomes will be anticorrelated if 
Alice and Bob measure along the same direction. In flat space the notion of 'same direction' is well-defined. 
However, in curved space, whether two directions are 'the same' or not is a matter of pure convention, since 
the direction obtained from parallel transporting a reference frame from Alice to Bob is path dependent. 
Thus, the phrase 'Alice and Bob measure along the same direction' does not have an unambiguous meaning 
in curved space. 

With no natural way to determine that two reference frames at separated points have the same orientation, 
we are left with having to keep track of the arbitrary local choice of reference frame at each point. The 
natural way to proceed is then to develop a formalism that will be reference frame covariant, with the 
empirical predictions (e.g. predicted probabilities) of the theory required to be manifestly reference frame 
invariant. The formalism obtained in this paper meets these two requirements. 

4.2 Tetrads and local Lorentz invariance 

The previous discussion was in terms of a curved space and a spatial reference frame consisting of three 
orthonormal spatial vectors. However, in this paper we consider curved spacctimes, and so we have to 
adjust the notion of a reference frame accordingly. We can do this by simply including the 4-velocity of the 
spatial reference frame as a fourth component t x of the reference frame. Thus, in relativity a reference frame 
(t x ,x Xl y x , z x ) at some point x consists of three orthonormal spacelike vectors and a timelike vector t x . 

Instead of using the cumbersome notation (t x , x x , y x , z x ) to represent a local reference frame at a point 
x we adopt the compact standard notation ej(x). Here / = 0, 1, 2, 3 labels the four orthonormal vectors of 
this reference frame such that ~ t, e^ ~ x, e% ~ y, and ~ z, and fi labels the four components of each 
vector with respect to the coordinates on the curved manifold. The object Cj(x) is called a tetrad field. This 
object represents a field of arbitrarily chosen orthonormal basis vectors for the tangent space for each point 
in the spacetime manifold Ai. This orthonormality is defined in spacetime by 

9fiu(x)e^(x)ej(x) = rju 

where g^ u is the spacetime metric tensor and rju is the local flat Minkowski metric. Furthermore, orthogo- 
nality implies that the determinant e = det(ej) of the tetrad as a matrix in (fx, I) must be non-zero. Thus 
there exists a unique inverse to the tetrad, denoted by e 1 , such that e^ej = T)j = Sj or e^ej" = g 1 ^ = 5". 
Making use of the inverse we obtain 

g^(x) = el{x)ei{x)n u . 

1 In this paper wc will implicitly always work in a topologically trivial open set. This allows us to ignore topological issues, 
e.g. the fact that not all manifolds will admit the existence of an everywhere non-singular field of reference frames. 
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Therefore, if we are given the inverse reference frame e^(ir) for all spacetime points x we can reconstruct the 
metric g llu {x). The tetrad e^{x) can therefore be regarded as a mathematical representation of the geometry. 

As stressed above, on a curved manifold the choice of reference frame at any specific point x is completely 
arbitrary. Consider then local, i.e. spacetime-dependent, transformations of the tetrad e^(x) — > e'i(x) — 
A/(x)e^(x) that preserve orthonormality; 

= g^(x)A I K (x)e^ K (x)A J L (x)el(x) = mL Aj K (x)Aj L (x). (4.1) 

The transformations A/ (x) are recognized as local Lorentz transformations and leave r\i j invariant. Given 
that the matrices A/ (x) are allowed to depend on so that different transformations can be performed at 
different points on the manifold, the reference frames associated with different points are therefore allowed 
to be changed in an uncorrelated manner. However for continuity reasons we will restrict Aj J (x) to local 
proper Lorentz transformations, i.e. members of SO + (l, 3). 

The inverse tetrad transforms as — > e'* — A I J e fl where A I K A J K — Sj. We now see that the 
gravitational field is invariant under these transformations: 

5^ = ^/je'/e'/ = ^jAVp A'l^ = vu^k^l^^ = VKLe^e^ = g^. (4.2) 

Therefore, all tetrads related by a local Lorentz transformation A 1 j(x) represent the same geometry 
g^ v . Thus, by switching from a metric representation to a tetrad representation we have introduced a new 
invariance: local Lorentz invariance. 

As stated earlier it will be useful to formulate qubits in curved spacetime in a reference frame covariant 
manner. To do so we need to be able to represent spacetime vectors with respect to the tetrads and not the 
coordinates. A spacetime vector V expressed in terms of the coordinates will carry the coordinate index V M . 
However, the vector could likewise be expressed in terms of the tetrad basis, in this case = V I ej where 
V 1 are the components of the vector in the tetrad basis given by V 1 = e^V^ 1 . We can therefore work with 
tensors represented either in the coordinate basis labelled by Greek indices /x, v, p, etc or in the tetrad basis 
where tensors are labelled with capital Roman indices /, J, K, etc. The indices are raised or lowered either 
with g pv or with i] IJ depending on the basis [^] We will switch between tetrad and coordinate indices freely 
throughout this paper. 



4.3 The connection 1-form 

In order to define a covariant derivative and parallel transport one needs a connection. When this connection 
is expressed in the coordinate basis, which is in general neither normalized nor orthogonal, this is referred to 
as the affine connection T p v . Alternatively if the connection is expressed in terms of the orthonormal tetrad 
basis it is called the connection one-form w„ j. To see this, consider the parallel transport of a vector 
along some path x ll (X) given by the equation 

DW _ dV» dx» 
DX ~ dA dA vp 

where A is some arbitrary parameter. The vector V 1 in the tetrad basis is expressed as V 1 = V 1 e^. We can 
now re-express the parallel transport equation in terms of the tetrad components V 1 : 

Di^V 1 ) _ d(e»V I ) dx» j 
' DX ~ dA + d\ Lv " ejV 



(dV 1 dx" 
\dX + ~dX 



Thus, if we define 



2 see the notation and conventions Section [] 
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the equation for the parallel transport of the tetrad components V 1 can be written as 



DV 1 _ dV 1 dx u j 



DX ~ dA ' dA^ J = °- 

The object ujJ j is called the connection 1-form or spin-1 connection and is merely the affine connection T^ p 
expressed in a local orthonormal frame e^{x). It is also called a Lie-algebra -valued 1-form since, when viewed 
as a matrix (uj^Yj, it is a 1-form in v of elements of the Lie algebra so(l, 3). The connection 1-form encodes 
the spacetime curvature but unlike the affine connection it transforms in a covariant way (as a covariant 
vector, or in a different language, as a 1-form) under coordinate transformations, due to it having a single 
coordinate index v. However, as can readily be checked from the definition, it transforms inhomogeneously 
under a change of tetrad e*(x) —> A 1 J {x)e'j 1 {x): 

<V J ~> u'Jj = A^A>/ L + AV„A/. (4.3) 

The inhomogeneous term A I K d fl A J K is present only when the rotations depend on the position coordinate 
x M and ensures that the parallel transport D ^ transforms properly as a contravariant vector under local 
Lorentz transformations. 



5 The qubit as the spin of a massive fermion 

A specific physical realization of a qubit is the spin of a massive fermion such as an electron. An electron 
can be thought of as a spin-^ gyroscope, where a rotation of 27r around some axis produces the original 
state but with a minus sign. Such an object is usually taken to be represented by a four-component Dirac 
field, which constitutes a reducible spin-| representation of the Lorentz group. However, given that we are 
after a qubit and therefore a two-dimensional object, we will work with a two-component Weyl spinor field 
4>a(x), with A = 1,2, which is the left-handed component of the Dirac field (see |A| . [^] The Weyl spinor itself 
constitutes a finite-dimensional faithful - and therefore non-unitary - representation of the Lorentz group 
|16j and one may therefore think that it could not mathematically represent a quantum state. As we shall 
see, unitarity is recovered by correctly identifying a suitable inner product. 

We will begin by considering the Dirac equation in curved spacetime minimally coupled to an electro- 
magnetic field. We rewrite this Dirac equation in second-order form (called the Van der Waerden equation) 
where the basic field is now a left-handed Weyl spinor <\>a- This equation is then studied in the WKB limit 
which separates the spin from the spatial degrees of freedom. We then localize this field along a classical 
trajectory to arrive at a transport equation for the spin of the fermion which forms the physical realization 
of the qubit. We find that this transport equation corresponds to the Fermi- Walker transport of the spin 
along a non-geodesic trajectory plus an additional precession of the fermion's spin due to the presence of 
local magnetic fields. We will see that from the WKB approximation a natural inner product for the two- 



dimensional vector space of Weyl spinors emerges. Furthermore, we will see in section 5.3.1 that in the rest 
frame of the qubit the standard notion of unitarity is regained. It is also in this frame where the transport 
equation is identical to the result obtained in |10) . 

5.1 The WKB approximation 

Before we begin our analysis of the Dirac equation in the WKB limit we refer the reader to [A] for notation 
and background material on spinors. This material is necessary for the relativistic treatment of massive 
fcrmions. 

5.1.1 The minimally coupled Dirac field in curved spacetime 

Fermions in flat spacetime are governed by the Dirac equation i^d^ — m^. Since we are dealing with 
curved spacetimes we must generalize the Dirac equation to include these situations. This is done as usual 
through minimal coupling by replacing the partial derivatives by covariant derivatives. The covariant deriva- 
tive of a Dirac spinor is defined by [UJ 

= (8^ - '-u^jS")* (5.1) 



3 We could work instead with the right-handed component, but this would yield the same results. 
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where S IJ = ^["f 1 ,J J ] are the spin-i generators of the Lorentz group and 7 7 are the Dirac 7-matrices which 
come with a tetrad rather than a tensor index. The gravitational field enters through the spin-1 connection 
Ufiij. We assume that the fermion is electrically charged and include an electromagnetic field Fjj by minimal 
coupling so that we can consider accelerated trajectories. The Dirac equation in curved spacetime minimally 
coupled to an external electromagnetic field is then given by 

i^D^ = to* (5.2) 

where we define the U(l) covariant derivative as D„ = V M — ieA^. 



5.1.2 The Van der Waerden equation: an equivalent second order formulation 

In order to proceed with the WKB approximation it is convenient to put the Dirac equation into a second- 
order form. This can be done by making use of the Weyl representation of the 7-matrices (see [A] for further 
details). In this representation the 7-matrices take on the form 

l_( a\ A , 



7 - 1 -I A' A q 

The Dirac equation then splits into two separate equations 

ia» A ' A D^ A =m X A ' (5.3a) 

\a^ AA ,D^x A = mc/)A (5.3b) 

with a^ A A = ej& IA A and o-^ A , A = e^a 1 ' A , A , and * = {4>AtX A )j where (f) A and \ A are left- and right- 
handed 2-spinors respectively. Solving for \ A m equation ( |5.3a[ ) and inserting the result into (5.3b I yields 
a second-order equation called the Van der Waerden equation |42j 

a il AA ,a vA ' B D tl D v 4> B + m 2 <t> A = 



which is equivalent to the Dirac equation (5.2). We can rewrite this equation in the following way 



= (rV^'^^fe+m 2 ^ 



g^D^D^ A - iL'Y V*v<b ~ [eS B F ^c + m^ A (5.4) 



where we have used that 2D^D V \ = [Z?^,!),,] and 2D^D U ^ = {Z^,!^}, and a^a^ = g^v. We identify 

B 
v [ivA 



Ffj, u = 2V[ M A„] as the electromagnetic tensor and ^■ tlv A B< t ) B '■= 2 Vj^ V^j <j> A as a spin-| curvature 2-form 



associated with the left-handed spin-| connection ^oj^ijL ij a b , where L^ v = e^ejL IJ = ^0^,0^ are the 
left-handed spin-| generators related to the Dirac four-component representation by S^ v = © B} iV . We 
have tacitly assumed here that the connection is torsion-free. Torsion can be included (at least in the case 
of vanishing electromagnetic field) and will slightly modify the way the spin of the qubit changes when 
transported along a trajectory. We refer the reader to [331 HH US] for further details on torsion. 



5.1.3 The basic ansatz 

The starting point of the WKB approximation is to write the left-handed two-spinor field 4>a as 

and study the Van der Waerden equation in the limit e — > 0, where e is a convenient expansion parameter. 
Physically this means that we are studying solutions for which the phase is varying much faster than the 
complex amplitude if a- In the high frequency limit e — > the fermion will not 'feel' the presence of a finite 
electromagnetic field. We are therefore going to assume that as the frequency increases the strength of the 
electromagnetic field also increases. We thus assume that the electromagnetic potential is given by -A^. e is 
to be thought of as a 'dummy' parameter whose only role is to identify the different orders in an expansion. 
Once the different orders have been identified the value of e in any equation can be set to 1. 
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5.1.4 The Van der Waerden equation in the WKB limit 

Rewriting the Van der Waerden equation in terms of the new variables ipA and 8, and collecting terms of 
similar order in -, yields 



gTVpVytpA ~ iL^/Vi^jftpc + ~ QkrVptpA + ^V M fc" + leF^L^/ <p B ) - ^k^tp A + m 2 tp A = (5.5) 



where we define the momentum/ wavevector as the gauge invariant quantity fc M = V ^9 — eA^. 

If we assume that both the typical scale I over which ip A varies and the curvature scale 1Z are large 
compared to the scale A over which the phase varies (which is parameterized by e), the first two terms of 
(5.5) can be neglected. In the WKB limit the mass term represents a large number and is therefore treated 



as a 1/e 2 term. The remaining equations are then 



2k»V^ A + (fA^ak" + \eFa V lf ,V A B VB = 



k^k^ - m 2 = 0. 



(5.6a) 
(5.6b) 



5.1.5 Derivation of the spin transport equation and conserved current 



The dispersion relation (5.6b I implies that k is timelike. Furthermore, by taking the covariant derivative of 



the dispersion relation and assuming vanishing torsion 



V„(fc"/c„-m 2 ) 



2{k"V ll k v + ek^F^) = 



we readily see that the integral curves of u^{x) = k^(x)/m, defined by -r 
force law 

D 2 x^ Ax v _ u n 



dx" 



, satisfy the classical Lorentz 



(5.7) 



where 



D 2 x» 
Dt 2 



^ r V i ,u' 1 and u^u^ = 1. Thus, the integral curves of k^ are classical particle trajectories. 



To see the implications of the first equation (5.6a) we contract it with k fl a ilA <fiA' and add the result to 



its conjugate. Simplifying this sum with the use of (5.7) and the identity ([46j Eqn (2.85) pl9]) 



KA 1 AjJ.l B _ 1 ( m KI- JA'B m KJ-IA'B 



rj a 



KIJ -LA'B\ 
r.CT 



yields 



where ip 2 = u^a^ A A <p A 'tp A . Eq.(|5.8[) can also be rewritten as 



V^W) =0 

which tells us that we have a conserved energy density j M = 



(5.8) 
(5.9) 



-w 



with g = de t ff Mi/ 



Secondly, (5.8) yields V M fc M = — (2fc M V M y;)/V and when this is inserted back into (5.6a) we obtain 



VkPV^A + ieF^L^/^B = 0. 



By making use of the integral curves x^(t) we obtain the ordinary differential equation 

D ^ A 4-i e TP T IJ B ih -0 



(5.10) 



where -jjf- — ^,u tl uj tl ijL IJ A B il>B is the spin-^ parallel transport. Equation (5.10) governs the evolution 

of the normalized spinor ip A = Va/v along integral curves. Below ipA will assume the role of the qubit 
quantum state. 



^ip 2 has dimension L 
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5.2 Qubits, localization and transport 

The aim of this paper is to obtain a formalism for localized qubits. However, the WKB approximation 
does not guarantee that the fermion is spatially localized, i.e. the envelope <p(x) need not have compact 
support in a small region of space. In addition, even if the envelope initially is well-localized there is nothing 
preventing it from distorting and spreading, and becoming delocalized. We therefore need to make additional 
assumptions beyond the WKB approximation to guarantee the initial and continued localization of the qubit. 
As pointed out in §3.2] by restricting ourselves to localized envelopes we avoid the particle number ambiguity 
and can interpret the Dirac field as a one-particle quantum wavefunction. 



5.2.1 Localization 

Before we begin let us be a bit more precise as to what it means for a qubit to be 'localized'. In order to avoid 
the particle number ambiguity we know that the wavepacket size C has to be much less than the curvature 
scale 1Z. We also know from quantum field theory that it is not possible to localize a massive fermion to 
within its Compton wavelength A com = h/mc using only positive frequency modes. Mathematically we 
should then have A com < C <C TZ where C is the packet length in the rest frame of the fermion. If A com ~ 7Z 
the formalism of this paper will not be empirically correct. 

How well-localized a wavepacket is, is determined by the support of the envelope. Strictly speaking we 
know from quantum field theory that a localized state will always have exponential tails which cannot be 
made to vanish using only positive frequency modes. However, the effects of such tails are small and for the 
purpose of this paper we will neglect them and assume that the wavepacket has compact support. 

The equation that governs the evolution of the envelope within the WKB approximation is the continuity 



equation (5.9) 



V„(uV(aO) = 0. 

If we assume that the divergence of the velocity field u M is zero, i.e. V^u^ 1 = 0, the continuity equation 
reduces to 

V^uVOr)) = u^V^tp 2 + ^Vpti" = u^V^ 2 = 0, 
or, using the integral curves of u M , 



dr 



= 0. 



Thus, the shape of the envelope in the qubit's rest frame remains unchanged during the evolution. However, 
because of the uncertainty principle [47] , if the wavepacket has finite spatial extent it cannot simultaneously 
have a sharp momentum, and therefore the divergence in velocity cannot be exactly zero. We can then relax 
the assumption, since the only thing that we need to guarantee is that the final wavepacket is not significantly 
distorted compared to the original one. Since V M u M measures the rate of change of the rest-frame volume 
y 137 EH] we should require that 

(V^) « — 

where (V^m^) is the typical value of IV^M^j, and T r the proper time along some path T assumed to have 
finite length. If we combine this assumption of negligible divergence with the assumption that the envelope 
is initially localized so that the wavepacket size is smaller than the curvature scale, we can approximately 
regard the envelope as being rigidly transported while neither distorting nor spreading during its evolution. 

To further suppress the spatial degrees of freedom we need also an assumption about the two-component 
spinor ipA(x). This variable could vary significantly within the localized support of the envelope <p{x). 
However, as we want to attach a single qubit quantum state to each point along a trajectory we need 
to assume that iPa(x) only varies along the trajectory and not spatially. More precisely, we assume that 
tpA(t, x) — ipA(t) when we use local Lorentz coordinates (t, x) adapted to the rest frame of the particle. This 
implies that the wavepacket takes on the form 
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This form is not preserved for all reference frames since in other local Lorentz coordinates tpA will have 
spatial dependence. Nevertheless, if the packet is sufficiently localized and if) A varies slowly the wave-packet 
will approximately be separable in spin and position for most choices of local Lorentz coordinates. With 
these additional assumptions we have effectively 'frozen out' the spatial degrees of freedom of the wavepacket. 
The spinor ipA can now be thought of not as a function of spacetime i/ja(x) satisfying a partial differential 
equation, but rather as a spin state ^Pa{ t ) defined on a classical trajectory Y satisfying an ordinary differential 



equation (5.11). We can therefore effectively characterize the fermion for each r by a position x^fr), a 4- 
velocity dx M /dr = u M (r), and a spin iPa(t). Once we have identified the spin as a quantum state this will 
provide the realization of a localized qubit. 

5.2.2 The physical interpretation of WKB equations 



As discussed in section 3.2 if we restrict ourselves to localized wavepackets we can interpret 4>a(x) as one- 
particle excitations of the quantum field. This allows us to interpret the conserved current j^/m = J~gip 2, u ii 
as the probability current of a single particle. In this way we can provide a physical interpretation of the 
classical two-component spinor field (\>a{x) as a quantum wavefunction of a single particle. 



Next, let us examine the transport equation (5.10). The electromagnetic tensor Fu that appears in 
the term ieFjjL IJ /m can be decomposed into a component parallel to the timelike 4- velocity u 1 and a 
spacelike component perpendicular to u 1 using a covariant spatial projector hj = 5j — v}uj . We can then 
rewrite F[jL IJ A B as (2uiu K Fkj + hj K hj" Fkl)L ij a B ■ The first term corresponds to the electric field as 
defined in the rest frame, uiu k Fkj- This will produce an acceleration u^W^ui = aj = —^u j Fjj of the 
fermion as described by the Lorentz force equation ( |5.7[ ). The second term is recognized as the magnetic 
field experienced by the particle, i.e. the magnetic field as defined in the rest frame of the particle, i?Jj st . 
We thus obtain the transport equation for ipA'i 

^ - iu iaj L IJ A B ^B + i^B}fL" A B ^B = 0. (5.11) 

This has a simple physical interpretation. The third term represents the magnetic precession which is 
induced by the torque that the magnetic field exerts on the spin. This takes the usual form ^^-FijL^ — 
^^Bij^e 13 k a k = — tj^B • a if we express it in a tetrad co-moving with the particle, i.e. = u M . 
The two first terms represent the spin-half version of the Fermi-Walker derivative: 

D FW j>A _ D^A . TlJB , r - 19 . 
— ^ — = miajL A 4>b- (5.12) 

The presence of a Fermi- Walker derivative can be understood directly from physical considerations. Heuris- 
tically we understand the electron as a spin-| object, i.e. loosely as a quantum gyroscope. The transport of 
the orientation of an ordinary classical gyro is not governed by the parallel transport equation but rather, 
it is governed by a Fermi- Walker transport equation. The Fermi- Walker equation arises when we want to 
move a gyroscope along some spacetime path without applying any external torque |48j . [^] We thus identify 



( 5.12 ) as describing torque- free transport of the electron, resulting in the usual Thomas precession of the spin 
|33j . Finally, the parallel transport term Di/ja/Dt encodes the influence of gravity on the qubit, governed 
by the spin-1 connection uiJj- 

5.2.3 A summary of the WKB limit 

Let us summarize the results from the previous section. 

• The full wavepacket is written as </>a(%) = ^A{x)^p{x)e l6 ^ x \ 

• The current j^/m = y/gip 2 ^ is a conserved probability density. 

• The phase 9 and the vector potential define a field of 4- velocities = ~(V \fi — eA^). 

• The integral curves of m m are timelike and satisfy the classical Lorentz equation ma^ = eu l 'F^ v . 



5 At first one might think that this is just what the parallel transport equation achieves. However, this is only true for 
geodesic motion (a/ = 0), where the Fermi- Walker and parallel transport equations agree. 
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• The two-component spinor ip A ( r ) defined along some integral curve of satisfies the transport equa- 
tion 

^ - \ Ul ajL IJ A B ^ B + i^-h I K h J L F KL L IJ A B il, B = (5.13) 
which dictates how the spin is influenced by the presence of an electromagnetic and gravitational field. 

5.3 The quantum Hilbert space 

The spinor ip A € W (where W is a two dimensional complex vector space) could potentially encode a 
two dimensional quantum state. However, given that i^a constitutes a faithful and therefore non-unitary 
representation of the Lorentz group this identification might seem problematic. This issue is resolved by 
identifying a velocity-dependent inner product on the space W . In doing so we are able to promote W to 
a Hilbert space and so regard tp A as a quantum state. Let us now show how the two-component spinor ipA 
can be taken as a representation of the quantum state for a qubit, and that it does indeed evolve unitarily. 

5.3.1 The quantum state and inner product 

Although the space of two-component spinors W is a two-dimensional complex vector space, it is not a 
Hilbert space as there is no positive definite sesquilinear inner product defined a priori. However, in the 
above analysis of the Dirac field in the WKB limit the object = uia IA emerged naturally. Note that 
this object is simply the inner product for the Dirac field in the WKB limit and has the appropriate index 



structure of an inner product for a spinor space (see A.3|. Thus we take the inner product between two 
spinors and ip A to be given by 

<^ 2 ) = l*HW A = u^'HWa (5-14) 

which in the rest frame u 1 = (1, 0, 0, 0) takes on the usual form uoa OA wa'^A = & A rA'^A- The connection 
between Dirac notation and spinor notation can therefore be identified as 

First note that the inner product ( |5.14[ ) is manifestly Lorentz invariant. This follows immediately from 
the fact that all indices have been contracted^] Secondly, I A A satisfies all the criteria for an inner product on 
a complex vector space W: Sesquilinearitjj^jis immediate, and the positive definiteness follows if u 1 is future 
causal and timelike, since the eigenvalues X± — u°(l±v) of I A A are strictly positive, where u° = (1 — v 2 )~z 
and v denotes the speed of the particle as measured in the tetrad frame. Thus, in the WKB limit, I A A can 
be taken to define an inner product on the spinor space W which therefore becomes a Hilbert space. The 
spinor i^a is then a member of a Hilbert space and thus it plays the role of a quantum state. A qubit is then 
characterized by its trajectory T and the quantum states tPa(j) attached to each point along the trajectory. 



In £2.1 we saw that we need a separate Hilbert space for each spacetime point x. However, the inner 
product is also velocity dependent, or equivalently momentum dependent. Thus, we must also regard states 
corresponding to qubits with different momenta as belonging to different Hilbert spaces. In particular, we 
cannot compare or add quantum states with different 4-momenta p\ ^ p2 even if the quantum states are 
associated with the same position in spacetime. Consequently the Hilbert space of the qubit is labelled not 
only with its spacetime position but also with its 4-momentum. We therefore denote the Hilbert space as 

5.3.2 Wigner rotation 

In order to establish a relation to the Wigner representations and Wigner rotations [T7] we first note that 
the basis 

U = (J) XA = (\ 



3 Lorentz invariance can be verified explicitly by making use of A 1 j(x)a JB B {x)hg, (x)hg(x) = a IA A [46| . 
Sesquilincarity is the property that the inner product is linear in its second argument and antilinear in its first. 
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in which the quantum state is expanded, ipA — ipi£,A + 'feXA; is an oblique basis and not orthonormal 
with respect to the inner product I A A , i.e. (£|x) 7^ and ^ 1 ^ (xlx)- One consequence of this is 
that the transport equation (5.12| appears non-unitary as it contains both terms that look Hermitianfe.g. 
L y = ^e l3 k a k ), and terms that look anti-Hermitian (e.g. L - 7 = — ^<7 J ). However, as will shall see in 98 the 
transport is unitary with respect to the inner product I A ■ 

The connection to the Wigner formalism is seen by re-expressing the quantum state in an orthonormal 
basis. This is given by 



Xa = A/xs 



where is the spin-^ representation corresponding to the Lorentz transformation defined by uj — A/ 5®. 
Orthonormality follows from the Lorentz invariance of a JB B and the fact that £a and xa are orthonormal 
with respect to the inner product 8 . For example, £ and x are orthogonal which can be seen by making 
use of the invariance of a IA A : 



(fix) 



-IA'B 



ic'^gXDUio* 

C ,/ A C A D \J <r0 -IA'B 
Z,C'XdA a , A b A jdj(T 



£,cXdS 



BC' 







and we can in a similar way demonstrate that = (xlx) = 1- The components (^1,^2) are defined by 
ipA = ipi£,A + 1P2XA and can now be understood as the components ipA = ^ 1 a ^b of the spinor in the 
particle's rest frame. 

Given that in the rest frame the basis (£,AtXa) is indeed orthonormal, it is instructive to also express 
the Fermi-Walker transport in such a basis. By doing so we will not only see that the evolution is indeed 
unitary, but in addition we will make contact with the transport equation identified by |10j in which the 
authors made use of infinite-dimensional representations and the Wigner rotations. 

Explicitly the spin-| Lorentz boost as defined above takes the form [46] 



A u - 

A A — 



7 + 1 
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a ~i b 



(5.15) 



where 6 l is the boost velocity, 7 = (1 — j3 2 ) 2 is its Lorentz factor, and the Pauli operators are given by 



T J B 



AA' 



-IA'B 



The corresponding spin-1 boost (acting on a contravariant vector) is 



A',= 



7ft 

r 7+1 



(5.16) 



where 



5ij0 l . Substituting ipA = A B ipg into the Fermi-Walker derivative (5.12| yields 



D FW ^ A dip. 



Dt 



A 

dr 



-u^w^ i.jL IJ A B ip B - iujdjL 



IJ B 

A 



Ipl 



= A, 



■Alps 
dr 



dAf - . 
dr 



ujaj ) L IJ A B A^ C = 0. 



The latter expression can be rearranged to give an evolution equation for the rest-frame spinor 



dtpA 
dr 



-A 



_ x sdA 



D 
B 



dr 



+ i I -u^u^ u + maj 



A-l «r/J Ca D 
A A ^ B A C 



IpD- 



One can then simplify this using the identities A B A D C L IJ 1 



D 



\I a .7 tKL 
A K A L A 







and A 



a b m P 9] 



to obtain spin-1 boosts to the terms involving L IJ . Using now explicit expressions of the spin-^ and spin-1 
boosts (5.15) and (5.16), one can cancel many of the terms to yield the result 



dr 2( 7 + iy l dr 



Qw M ij + 7^ ojPi + uP l Pj \ L%3 a$b- (5.17) 



3 This can be shown using the Lorentz invariance of a IA A and the definition of L IJ a in terms of a 1 : see [A. 2 
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This is the transport equation for the quantum state i\>a expressed in terms of the rest- frame spinor ip^. First 
we note that the transport is unitary with respect to the standard inner product S A A as it only contains 
terms proportional to = \e lJ k a k . It is however not manifestly Lorentz invariant. Secondly, it is also 
equivalent to the transport equation derived by |10j who used the infinite-dimensional Wigner representations 
[17] . We have thus re-derived their result using the Dirac equation in the WKB limit. In addition, we have 
done so while avoiding the use of momentum eigenstates \p, a) , which are strictly speaking not well-defined 
in a curved spacetime as no translational invariance is present. 

Notice that there is no term proportional to the identity 5P which would correspond to an accumulation 
of global phase. In fact, global phase is missing in [TO]. On the other hand, as we shall see in ^J7j these 
phases are automatically included in the WKB approach adopted in this paper. 

Although the unitarity of the transport becomes manifest when written in terms of the rest-frame spinor 
it is not necessary to work with equation (5.17). Once we generalize the notion of unitarity in f|8]we will 
see that we can treat the evolution of the quantum state in terms of the manifestly Lorentz covariant 
Fermi-Walker transport (5.12). 



6 The qubit as the polarization of a photon 

Another specific physical realization of a qubit is the polarization of a single photon. This is an important 
example since it lends itself easily to physical applications. We obtain this realization via the WKB limit 
of Maxwell's equations in curved spacetime [3H1 EH]- The polarization of a photon is described by a unit 
spacelike 4- vector ip^ called the polarization vector [JHlEniSS] ■ Restricting ourselves to localized wavepackets 
we obtain the description of a photon with definite 4-momcntum/wavevector k^ 1 and polarization vector 
which is parallel transported along a null geodesic x M (A). We will see that in fact Vv contains only two gauge 
invariant degrees of freedom and thus can be taken to encode the quantum state of a photonic qubit. 

Although we consider only geodesic trajectories in this paper it is possible to consider non-geodesic tra- 
jectories. We refer the reader to [B] for a discussion of approaches to this problem. A physically motivated 
way to obtain non-geodesic trajectories would be to introduce a medium in Maxwell's equations through 
which the photon propagates. Nevertheless, even without explicitly including a medium, it is easy to in- 
clude optical elements such as mirrors, prisms, and other unitary transformations as long as their effect on 
polarization can be considered separately to the effect of transport through curved spacetime. 



6.1 Parallel transport from the WKB approximation 

In this section we shall see that the parallel transport equation for the polarization vector emerges directly 
from the WKB approximation (48J [49] . Gauge invariance and gauge fixing in the WKB approach are 
important for properly isolating the quantum state and we have therefore paid attention to this issue. 



6.1.1 The basic ansatz 

The WKB approximation for photons follows a procedure similar to that for the Dirac field. First we write 
the vector potential as 

A»=M<P^ ,e \- (6-1) 

As in the case for the Dirac field, the WKB limit is where the phase 9 is oscillating rapidly compared to 
the slowly varying complex amplitude tp^. As before, this is expressed through the expansion parameter e. 
Maxwell's equations can then be studied in the limit e 0. Although we omit taking the real part of <p M e ie / e 
in this section it is implicitly understood that this is done. 



6.1.2 Gauge transformations in the WKB limit 

Let us now study the U(l) gauge transformations in terms of the new variables 9 and ip^. It is clear that 
not all gauge transformations — > A^ + V^A will preserve the basic form A^ = ip^e 10 ^ 5 . We therefore 
consider gauge transformations of the form A = Qe l8 / e where £ is a slowly varying function. This class of 



gauge transformations can be written in the polar form of (6.1 1 as 



A^ -> A^ + V M A = A^ + V M (Ce ie / £ ) = + V P C + e : 
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and so ip^ -> ^ + V^C + ;^C- 

In the limit e — > note that <p^ docs not behave properly under the gauge transformations of the type 
that we are considering since the second term blows up. This has no physical significance and is just an 



artefact of describing the vector potential as being of the specific form (6.1). Such a gauge transformation 



leaves the physics unchanged but will no longer preserve the form of the solution (6.1) where we have a 



slowly varying envelope and rapid phase. Because of this it is necessary to further restrict the space of gauge 
transformations to "small" gauge transformations £ = — ie£. In that limit we then have 



Pv^Kfij.- ieV M £ + k^ 



(6.2) 



and so p>^ — > p>^ + + 0(e). However, as we shall see below, in order to maintain gauge invariance of the 
equations in all orders of e it is important to keep both orders of e in the gauge transformation (6.2). 



6.1.3 The gauge condition 

In the literature we find two suggestions for imposing a gauge. For example, in |48| the Lorenz gauge is 
used, V M A'' = (V '^(p 1 * + i^p^Q'S/t — o, and in the gauge fc'V^ = is imposed so that the complex 
amplitude ip^ is always orthogonal to the wavevector k^. However, for our purposes neither of these gauge 
conditions turns out to be suitable. Rather we will work in a gauge where k^ and ip^ are orthogonal up to 
first-order terms in e, i.e. 

p^ = ea(x) 
where a is taken to be some arbitrary function of x^. 



6.1.4 Maxwell's equations in the WKB limit 

Let us now turn to Maxwell's equations in vacuum: 



0. 



(6.3) 



The equations V^-F^j,] = are mere identities when we work with a vector potential rather than the 

al6l / £ into (6.3) we obtain 



gauge invariant = V^.A„ — V„A^ t . If we substitute the ansatz A p 



1 



D<p v - VV^ + ~(2fc"V^ + VvV^ ~ tV M / - V^fc")) - ^(fcV - k^X) = 0- ( 6 - 4 ) 



Gauge invariance can be a bit subtle in this context so let us make a few remarks. Eq.(6.4) is of course 
invariant under gauge transformations p>^ — >• tp^ — ieV + k^ as this is nothing but Maxwell's equations 
(6.3) rewritten in different variables. However, note that the terms of zeroth, first, and second order (in 1/e) 



ofEq. (6.4) 



2fc^V Al </?„ + ^V P F 
k 2 p v - k u ip^ 



(6.5a) 
(6.5b) 
(6.5c) 



are not separately gauge invariant. This is so because the gauge transformation i/3„ — >• i/3„ — ieV M £ + 



contains terms of different orders in e. Thus, after a gauge transformation of the second-order term (6.5c ) we 



end up with first-order terms in e, which then belong to (6.5b). Similarly first-order terms in e in (6.5b) end 



up in (6.5a). It is then easy to verify that the entire equation (6.4) is gauge invariant although the separate 



terms in (6.5) are not. 



6.1.5 Equations of motions in the gauge p^k* 1 = ea 



Imposing the gauge condition k^ip^ — ea on (6.4) yields the equation 



U<p v - W7 v <p M - V„a + -(2A^V^ + p v V ^ - fc„(V ^ - a)) - -^k 2 p u 



(6.6) 
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We now demand that the solutions for (p^ be independent of e in the limit when e is small. Physically 
this means that for high frequencies the form of the solutions should be independent of the frequency 
(parameterized by e). Consequently, each separate order of - in the expansion must be zero. The equations 
corresponding to the first and second orders then read 



k^k^ = 



(6.7a) 
(6.7b) 



for ip v ^ 0. The zeroth-ordcr equation is to be thought of as 'small' in comparison to the higher order terms 



in 1/e and is therefore ignored and not imposed as an equation of motion. The second equation (6.7b) is 
trivially gauge invariant since k^ does not transform. The first equation is only gauge invariant up to first- 
order terms in e. This can be seen by letting a transform as a — > a + fc M under a gauge transformation, 
making use of (6.5c), and the fact that k^ satisfies the geodesic equation as shown in (6.8 



6.1.6 The derivation of parallel transport and conserved currents 



Equation (6.7b) tells us that the wavevector fc M is a null vector, and that its integral curves x M (A) defined 
by dx^/dX oc fc M lie on a light cone. Taking the derivative of Eq.(6.7b) yields 



V v {k^k^) = 2k< t V v k ll = 2*^^(9 = 2fc"V ([1 V I/ = 2fc"V p fe I , = 



(6.8) 



which tells us that the integral curves are null geodesies]^] These are expected since we have considered 
Maxwell's equations in vacuum. Non-geodesic trajectories can be obtained by introducing a medium through 
which the photon propagates. Sec [B] for a discussion. 



Contracting equation (6.7a) with ip v and adding to it its complex conjugate yields the continuity equation 

2(p v k^^ v + 2^fc M V M ^„ + 2Cp v ip u V ^ = -2V A1 (<^ 2 fc' i ) = (6.9) 

where ip 2 = —g^'Pfj.fv Note that f 2 is gauge invariant up to first-order terms in e, i.e. tp 2 — > (p 2 + 0(e), 
and therefore also j M = ^fcjip 2 ^ is gauge invariant to first order. This means that is a conserved current 
in the WKB limit. Since j° has the units of a probability densitj{^] we can interpret j M as a conserved 
probability density current. 



We can also deduce that V^/c^ = — ^k^V^ip and if we insert this in equation (6.7a) and define the 



polarization vector i\> v through tp v = (pipv, we obtain 



2k»V^ u + ipvV^ - k y {V^ -a) = 2/c^V^ - ^-A"V M v - ^ 



a) 



which implies that 



2ipk»V^ v - k v {V^ - a) = 



2 r 



However, since a is arbitrary the whole right-hand side is arbitrary and we can write 



(6.10) 



The right-hand side is proportional to the wavevector k v and represents an arbitrary infinitesimal gauge 
transformation of ip^. Let us now introduce the integral curves of — k^/E given by dx^/dX — it' 1 where 
E is an arbitrary constant with dimensions of energy. We can then write equation (6.10) as 



DX 



(6.11) 



Thus the transport of the polarization vector is given by the parallel transport along the null geodesic 
integral curves of u^, with an arbitrary infinitesimal gauge transformation at each instant. 

'We have assumed in l |6.5a| l that the spacetime torsion is zero. A non-zero torsion field could possibly influence the 
polarization (see |45|). 

10 We recall that A„ has dimensions L" 1 in natural units with e = 1. 
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6.2 Localization of the qubit 

As in the fermion case, the WKB approximation is not enough to guarantee either that the wavepacket is 
localized or that it stays localized under evolution, and again it is not possible to achieve strict localization. 
Indeed it can be proved that a photon must have non- vanishing sub-exponential tails |31| 132) . As in the 
case of fermions, §5.2| we are going to ignore these small tails and treat the wavepacket as effectively having 
compact support within some small region much smaller than the typical curvature scale. 



The continuity equation (6.9) dictates the evolution of the envelope <p(x). Divided by the energy as 
measured in some arbitrary frame it becomes V At (</3 2 M A ') = 0. Again we see that the assumption V M u M = 
simplifies this equation. However, the interpretation of V^it^ 1 is a bit different. Instead of quantifying how 
much a spatial volume element is changing (as in j|5.2[), it quantifies how much an area element, transverse 



to in some arbitrary reference frame, changes |51|: 

rv u 1 dA 
V ^ = 2dA 

where A is an affine parameter defined by dx M /dA = u^. In this case we require that (V^u^) <C 1/Ar, 
where Ar is the affine length of the trajectory T. Thus, it gives us a measure of the transverse distortion 
of a wavepacket. For photons there can be no longitudinal distortion since all components, regardless of 
frequency, travel with the speed of light. Initial localization and the assumption that V^u^ w therefore 
guarantee that the wave-packet is rigidly transported along the trajectory. 

Once we assume that the polarization vector ipi does not vary spatially within the wavepacket we can 
effectively describe the system as a polarization vector ip^X) f° r eacn A G T. Having effectively suppressed 
the spatial degrees of freedom of the wavepacket, the polarization tp^ can thus be thought of as a function 



defined on a classical trajectory T, satisfying an ordinary differential equation (6.11). A photonic qubit can 
then be characterized by a position a;'' (A), a wavevector fc M (A), and a spacelike complex- valued polarization 
vector ip^X). 

6.3 A summary of WKB limit 

To summarize, the WKB approximation yields the following results and equations: 

• The integral curves ic M (A) of are null geodesies 

• The vector j M = y/g(f 2 k^ is a conserved probability density current. 

• The polarization vector tp^ satisfies tp^u^ — and transforms as tp^ — ► + vu v under gauge trans- 
formation up to first-order terms in e. 



• The transport of tpn is governed by (6.11 ) which is simply the parallel transport along integral curves 
of u M modulo gauge transformations. 

We have now established a formalism for the quantum state of a localized qubit which is invariant under 
ipfi —} tpfj, + vu^ and V-V uAl — up to first-order terms in e. We shall from this point on neglect the small 
terms of order e. 



6.4 The quantum state 

We now show that the polarization 4-vector has only two complex degrees of freedom and in fact it can be 
taken to encode a two-dimensional quantum state. We do this first with a tetrad adapted to the velocity 
of the photon for simplicity and then with a general tetrad. It is convenient and more transparent to work 
with tetrad indices instead of the ordinary tensor indices and we shall do so here. 

6.4.1 Identification of the quantum state with an adapted tetrad 

Recall from the previous section that we partially fixed the gauge to ujip 1 — 0. The remaining gauge 
transformations are of the form ip 1 ip 1 + vu 1 . Indeed, if ujip 1 = we also have that ujiijj 1 + vu 1 ) = for 
all complex- valued functions v, since u 1 is null. 
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To illustrate in more detail what effect this gauge transformation has on the polarization vector we adapt 
the tetrad reference frame ej to the direction of the photon so that u M oc eg + eg. Notice that there are 
several choices of tetrads that put the photon 4- velocity into this standard form. The two-parameter family 
of transformations relating these different tetrad choices are (1) spatial rotations around the z-axis and (2) 
boosts along the z-axis. 

With a suitable parameterization of the photon trajectory such that e^da^/dA) = 1 we can eliminate 
the proportionality factor and we have = e^ + eg . In tetrad components u 1 = (1, 0, 0, 1) and we see that 
the tetrad z-component eg is aligned with the photon's 3-velocity. Since = uj^p 1 — ip° — ip 3 it follows that 
ip° = ip 3 = v and the polarization vector can be written as 



iP 1 

w 



It is clear that a gauge transformation 



w 



vu 



(v + v\ 

ip 1 

\v + vj 



leaves the two middle components unchanged and changes only the zeroth and third components. The two 
complex components ip 1 an d ip 2 , which form the Jones vector |52j . therefore represent gauge invariant true 
degrees of freedom of the polarization vector whereas the zeroth and third components represent pure gauge. 

We can now identify the quantum state as the two gauge invariant middle components ip l and ip 2 , where 
ip 1 is the horizontal and ip 2 the vertical component of the quantum state in the linear polarization basis: 



1 



V o J 



|2) 





1 

V o / 



or simply | ^4} ~ 8 A with A — 1,2. The quantum state is then 



Note we have deliberately used a notation similar to that used for representing spinors; however, ip A should 
not be confused with an SL(2,C) spinor. In order to distinguish ip A from ip 1 we will refer to the former as 
the Jones vector and the latter as the polarization vector. 

6.4.2 Identification of the quantum state with a non-adapted tetrad 

In the above discussion we have used an adapted tetrad in order to identify the quantum state. We can 
write a map for this adaption explicitly, which will provide a generic non-adapted formalism. To adapt one 
simply introduces a rotation which takes the 4-velocity u 1 to the standard form [17] 



w 



which results in the tetrad being aligned with the photon's 3-velocity, as illustrated in figure 6.1 
Such a rotation is explicitly given by 



R z j(u) 



°0°J 



«3 
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Rei 




Figure 6.1: The rotation R adapts the spacelike vectors of the tetrad ei so that the z-axis e*3 is 
aligned to the 3-velocity of the photon u. A polarization vector is then in the plane spanned by 
Re\ and Re^. 



where u^/u — u^e^/u = — cos 9 is the angle between the direction of the photon and the z-component of 
the tetrad, r 1 = £j 30 u J is the spatial axis of rotation with f 1 = r 1 /\r\, and P I J = Sj — e\e®j is the projector 



onto the spacelike hypersurface orthogonal to the tetrad time axis (see Fig 6.1 ). 

It is important to stress that there are several other possible choices for this spatial rotation corresponding 
to different conventions for the linear polarization basis. Furthermore, the rotation matrix above becomes 
undefined for 8 — tt which is unavoidable for topological reasons. 

The rotation R J j induces a linear polarization basis SfB/j. We can now extract the components of the 
quantum state expressed in this basis as 

^ A =fjV, with ff = 5fB*j. (6.12) 

It is clear that the specific linear polarization basis used here depends on how we have adapted the tetrad 
to the velocity of the photon. However, regardless of what convention one chooses, the quantum state ip A is 
gauge invariant. Alternatively we could think of the quantum state directly in terms of an equivalence class 
of polarization vectors t/r ~ ip 1 + vuI orthogonal to photon velocity u 1 . The advantage of this approach 
is that once one has developed a gauge invariant formalism one need not work with the cumbersome two 
component Jones vector, but instead can work solely with the gauge covariant polarization vector tp 1 . This 
will be addressed below and in Section [8T2l 

ff from (6.12) turns out to provide a 'diad' frame: The two vectors fj and fj span the two-dimensional 



space orthogonal to both the photon's 4- velocity uj and the time component of a tetrad e\. If we let f\ 
be the inverse of ff we have that f^ff = Sj and fj^Uj = f^e^ — 0. In fact, if we define wj to be a null 
vector defined by e\ = \{ui + Wi) [ST], the vectors uj, wj, fj and fj span the full tangent space. This 



decomposition will be useful when identifying unitary operations in Section 8.2 



6.5 The inner product 

We must identify an inner product on the complex vector space for polarization so that it can be promoted 
to a Hilbert space. In the analysis of the WKB limit we found that j 1 = — % /gk I ip 2 rjjKil) J '4> K corresponded 
to a conserved 4-current which was physically interpreted as a conserved probability density current. A 
natural inner product between two polarization 4-vectors ip 1 and 4> J is then given by 

- WV 7 . (6.13) 

This form is clearly sesquilinear and positive definite for spacelike polarization vectors^] Unlike the case for 
fermions, the inner product r\u is not explicitly dependent on the photon 4-velocity. However, if we consider 
the gauge transformation ijj 1 — > tfj 1 + viu{ and (j) 1 — > (j) 1 + viu\ it is clear that unless u[ — uf, i.e. k[ cx 



the inner product (6.13) is not gauge invariant. We conclude that two polarization vectors corresponding to 
two photons with non-parallel null velocities do not lie in the same Hilbert space. Furthermore, in order to 
be able to coherently add two polarization states it is also necessary to have k[ = i.e. the two photons 
must have the same frequency. Under such conditions the inner product is both Lorentz invariant and gauge 



invariant. With the inner product (6.13) the complex vector space of polarization vectors is promoted to a 
Hilbert space which is notably labelled again with both position and 4-momentum p 1 = Hk 1 . 



11 There is no primed index for conjugate terms because the vector representation of the Lorentz group is real. 
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The above inner product (6.131 reduces to the standard inner product for a two-dimensional Hilbert 



space. This is best seen through the use of an adapted tetrad. In an adapted frame the inner product of 
ip 1 = (v 1 ip 1 ,ip 2 , v) with some other polarization vector (p 1 — (/i, (p 1 , (p 2 , /i) is given by 



-r)u(j> tp 



-flV + (pilp 1 + folp 2 + [IV = Cpilp 1 + 4> 2 ip 2 = (4>\ip). 



Thus, the standard inner product (<p\4P) — (piip 1 + <p2' l P 2 is simply given by (<p\ip) = —r/jj(p I ip J , where we 
associate 



ip and 



-ViJ0 J 



We can now work directly with the polarization 4-vector tp 1 which transforms in a manifestly Lorentz 
covariant and gauge covariant manner. 



6.6 The relation to the Wigner formalism 

The Wigner rotation W A B (k, A) on the quantum state represented by the Jones vector which results from 
the transport of the polarization vector can be identified in the same way as was done in |5.3.2| for fermions. 
Specifically, this is achieved by determining the evolution of the Jones vector ip A that is induced by the 
transport of the quantum state represented by the polarization 4-vector ip 1 . Substituting ip 1 = f B ip B m t ne 
transport equation (6.11 ), + u^uj 1 jip J — piu 1 , and multiplying by ff, we obtain 



dip A 



-V j 



f J B + h 



dA 



(6.14) 



The last term is zero, as ff is the diad frame defined to be orthogonal to u 1 . If we first consider (6.14) 



in an adapted tetrad as in ( 6.4.1 



right-hand side can be simplmcdto 



we see that the derivative -4f- vanishes and the remaining term on the 



dip A 
~d~X 



A j.B 



(6.15) 



where we have made use of the antisymmetry of the spin-1 connection in order to introduce the antisymmetric 
Pauli Y matrix cr^^J^] Eq.(6.15| is then clearly unitary and helicity preserving as it is proportional to 
<j v B in the linear polarization basis. We can now readily identify the infinitesimal Wigner rotation as 
iit /J ''io ^nu 1 ' B , where the rotation angle is u^uj^. In a non-adapted tetrad frame the map j l A 

f-l J 



W A B 



can be seen to put (|6.14| in the form (|6.15|) with a modified spin-1 connection ui' u i 1 . The Wigner rotation 



for non-adapted tetrads is then 



W A B = iu' t (i?/a M i? 2 J + iZ/w 



R 2 j)<r yA B 



(6.16) 



A general Wigner rotation is understood as the composition of maps W A B = f B (Au)A/f A (u). It is therefore 
no surprise that the transport of the polarization vector induces a Wigner rotation: The action of the 
gravitational field alon g a trajectory is simply a sequence of infinitesimal Lorentz transformations which are 
given by u^utj j ({ 2.2 ). The transport of the Jones vector is therefore described by a sequence of infinitesimal 
Wigner rotations given by (6.16). Notice that the Wigner rotation takes on a form which is not manifestly 



Lorentz covariant. This is because the Wigner rotation describes a spatial rotation. This should be contrasted 
with the manifestly Lorentz covariant representation in terms of parallel transported polarization vectors. 
Furthermore, after we have developed a measurement formalism in section |8.2.2| it will become clear that 
there is no need to work with the cumbersome Wigner rotations. 



7 Phases and interferometry 



So far we have determined the transport of the quantum state of a single qubit along one spacetime trajectory. 



If we inspect the transport equations (5.11) and (6.10) we see that neither one contains a term proportional 



2 Note again that this should not be confused with the cr-matrices encountered when working with spinors. 
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to the identity (6^ for fermions and Sj for photons). Such a term would lead to an overall accumulation 
of global phase e l6 ipA or e^ip 1 . This leads one to suspect that not all of the possible contributions to the 
global phase have been taken into account in these transport equations. Indeed this is the case, as can be 
seen immediately by considering the full wavepacket in the WKB approximation 

where <p(x) is a real-valued envelope and a = 1,2 for fermions or a = 0,1,2,3 for photons p^] Clearly there 
is an additional phase 6(x) which is not included in ip a {x). 

Since global phase is unobservable this is of course of no concern if we restrict ourselves to a qubit moving 
along a single trajectory. However, quantum mechanics allows for more exotic experiments where a single 
qubit is split up into a spatial superposition, simultaneously transported along multiple distinct paths, and 
recombined so as to produce quantum interference phenomena. Here it becomes necessary to keep track of 
the phase difference between the components of the spatial superposition in order to be able to predict the 
measurement probabilities at the detectors. 

In this section we will extend the formalism in this paper to include gravitationally induced phase 
difference in experiments involving path superpositions. The formalism will be derived from equations of 
the WKB approximation together with the assumptions of localization. With these assumptions, the details 
of the spatial profile of the qubit become irrelevant, and we can satisfactorily describe the experiment solely 
in terms of a phase difference A9 between two quantum states. This phase difference will depend on the 
spacetime geometry g^ v and the trajectories along which the qubit is simultaneously transported. We show 
how the various sources for the phase difference can be understood from a wave-geometric picture. We 
lastly apply the formalism of this paper to gravitational neutron interferometry [271 1261 125] and obtain an 
exact general relativistic expression for the phase difference which in various limits reproduces the results in 
|23[ I24[ I29j in which higher order corrections to the non-relativistic result were proposed. 



7.1 Spacetime Mach— Zehnder interferometry 

We consider, as a concrete example of an interference experiment, standard Mach Zehnder interferometry. 
As usual, there is a qubit incident on a beam splitter (e.g. a half-silvered mirror) which creates a spatial 
superposition of the qubit. The two components of the spatially superposed state (each assumed to be 
spatially well-localized) are then transported along two different paths and later made to interfere using 
another beam splitter. This produces two output rays each incident on a particle detector, as illustrated in 



figure 7.1 



Let us now focus our attention on a small region D\ situated on the right output arm. There are then 
two classical paths Y\ and r 2 which arrive at D±, as illustrated in Fig |7.1[ Note that in order for us to derive 
a formalism in terms of quantum states, the wavevectors and fc 2 of I^ and T 2 must be approximately 
equal in this region, i.e. k± = k% = k^. This is because the Hilbert space of a quantum state is labelled with 



momentum, as explained in £5.3.1 and ^6.5 



In general the times of arrival of the two paths Ti and F 2 at the second beam splitter BS2 will differ for 
the two paths. As we shall see this contributes to the total phase difference between the two packets. 

Let x be some suitable local Lorentz coordinate system in region D\. The wavepacket in region D\ is 
then given by the superposition 

a¥ a 1 Hx) + b¥ a 2 \x) (7.1) 

where ^i? \x) and ^^\x) are the packets propagated along ri and T 2 respectively, a and b are determined 
from the reflection and transmission coefficients of the various beam splitters in the experiment. In the case 



of 50-50 beam splitters, a = b = -4= in region D\ (see Fig 7.1 ). We will ignore any overall global phase factor 
resulting from reflections. 



7.2 The phase difference from the WKB approximation 

In order to make empirical predictions in a Mach-Zehnder type interference experiment we must determine 
explicitly the forms of ¥„\x) and ^a{x) in (7.1) in the detector region D±. By making use of the field 



13 In the case of a scalar particle (and thus not a qubit) there are still gravitational phases, and here the index a can just be 
removed. 
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Figure 7.1: Spacetime figure of a Mach-Zehnder type interferometer, illustrating a single qubit 
subjected to a beam splitter BS\ resulting in a superposition of the qubit travelling along two 
distinct spacetime paths Fi and IV In some future spacetime region containing the beam splitter 
BS2 the components of the spatial superposition are assumed to recombine to produce possible 
interference phenomena in the detector regions Di and D2. Fhese regions contain two trajectories, 
indicating that the times of arrival at the second beam splitter BS2 are not in general the same. 
The variables a?i,2 are arbitrary spacetime points in the region D\ along trajectories Ti and T2 
and are useful for calculating the total phase difference. The red and blue strips represent the 
spatial extents of the wavepackets along Ti and F2 respectively. These correspond to the length 
C of the wavepacket as measured along the line joining the points x\ and xi- 



equations in the WKB limit and the localization assumptions we will see that v£\x) and ^ {x) will 
differ by a phase accumulated along the trajectory and a rigid translation/displacement, resulting in an 
overall phase difference. The derivations differ in the cases of fermions and photons and we will treat them 
separately. 

7.2.1 Fermions 

In the small region D\ the wavepacket in the WKB approximation is given by 

a^\x) + bcj ) { 2 ) (x) = a Vl {x)^\x)e w ^ + b^x)^ (x)e ie ^ (7.2) 
where x is some local Lorentz coordinate system, and a an d b are real-valued coefficients. The functions 



<t>i{x), ^(x) and 6i{x) (i = 1,2) are defined in section 5.1 We are now going to successively make use of 
the equations of the WKB approximation and the localization assumptions to simplify the expression (7.2) 
and thereby extract the relative phase difference between the two components in the superposition. 

First we use the fact that under the mathematical assumptions detailed in §5.2. 1| the envelope will 
be transported rigidly and will not distort. Therefore, <fi(x) and ip2 (x) will differ at most up to a rigid 
translation and rotation. We assume that the packet is 'cigar shaped' and is always oriented in the direction 
of motion. The final envelopes will then differ at most up to a translation and we can write (fi(x) = ip(x — Xi) 
for some suitable function ip(x) and an arbitrary choice of spacetime points Xi,x% G D\ situated on the 



trajectories Fi,r 2 respectively (see figure 7.1| 
If we now assume that 



(p{x) 



for all points x € X>i, the difference in the envelopes (fii(x) and f2(x) is negligible. The translational 
difference in the envelopes can then be neglected and factored out: 

a4> [ l\x) + b(b { 2\x) w p(a;) (aV^riK 61 ^ + b^\r 2 )e w ^ (7.3) 
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where ip^(Ti) are determined by integrating the transport equation (5.11) to the points x\, x%, respectively. 
Thus for the purpose of interfcromctry the details of the envelope become irrelevant and can be ignored. 



We now focus on the phase 9\ and 9 2 . As we pointed out in section 5.3.1 in order to coherently add two 
quantum states it is necessary to assume that the wavevectors of the packets are the same, = k 2 — k^. 
Therefore in the region Di we have from the WKB approximation that the phases 0\(x) and 6*2(2:) both 
satisfy the equation 



k,, + eA,, 



(7.4) 



Within the small region D\ we regard k^{x) and A^(x) as constant and so the partial differential equation 
(7.4) has the solution 9i{x) = (fc M + eA,j){x^ — x^) + 0i(xi) where Oi(xi) are two integration constants 
corresponding to the value of 9i{x) at the points xi. These integration constants can be determined by 
integrating (7.4) along the trajectories r 12 to the positions x\^ 2 respectively, i.e. 



Oi{xi) 



(7.5) 



Using the above we can rewrite (7.3) as 



where 9$ is some arbitrary global phase just before the wavepacket was split up by the first beam splitter. 

(7.6) 



a^\x) + b<j>f(x) « cp(x)e w ^ (aV^^i) + b^(T 2 )e iA6 ) 



where 



Ad = (fc„ + eA^xt - a#) + (9 2 (x 2 ) - 9 1 (x 1 )). 



(7.7) 



It is important to note that this phase difference is independent of x^ G D± as all dependence on x has been 



factored out in (7.6). Furthermore, we note that the choice of x\ and x 2 is arbitrary and the phase difference 
A9 is also independent of this choice. To see this, consider a different choice of positions, x\ = x\ + 8x1 
and x' 2 — X2 + Sx 2 on Ti and IV This results in a change in the integration constants (7.5) of 9i(xi) — > 
9i{xi) + (kn + eA^Sxj which exactly cancels the change in the term (fc M + eA fJ _)(x^ — x%) in ( 7.7 ). Therefore 



A9 is independent of the arbitrary positions x% and x 2 . 

Note that A9 is not the phase difference determined empirically in a Mach-Zehnder type interference 
experiment. This is because the transported quantum states ^>^(ri) and ip^'(T^) can contain an additional 
phase difference induced from their specific evolutions on the Bloch sphere. This transport induced phase 
difference can be determined from 1531 



AA0-J 



KVW(Ti)|V( a) (r a )>r 



(7.8 



The region D\ is assumed to be small enough that ijj^ and tp^ do not vary significantly with changes 
in xi and x 2 . Thus, the total phase difference A0Tot, which is the quantity that we actually measure in a 
Mach-Zehnder experiment, is then given by 

A6*Tot = A9 + A6*Trans- 

This total phase difference A6*Tot can be determined completely from the trajectories Ti and T 2 and the 



spacetime geometry using the transport equation (5.11). In particular, the phase difference measured 
by some detector in D\ is independent of the motion of that detector. 

Lastly, if we restrict ourselves to measurements that do not probe the spatial profile we can neglect the 



factor ip(x)e Wl ^ in (7.6). All contributions to the phase difference are then contained in ip^ and A9, and 



so at Di we are left with the two-dimensional quantum state 



>rccomb = a ^A } ( T l) 



b^ ] (T 2 )e 



1A8 



Therefore the assumptions that led to ( 7.6 1 established a formalism for determining the resulting qubit 
quantum state in region D\ of a Mach-Zehnder type interferometer. 
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7.2.2 Photons 

The derivation of the phase difference for photons follows essentially the same path as that for fermions. 
The starting point is to consider the wavepacket in the small region D\ 

aA\ x) {x) + bA\ 2) {x) = aip x {x)il>\ 1) {x)e ie ^ + b^ 2 (x)^ (2) {x)e l9 ^ (7.9) 

where x is some local Lorentz coordinate system, and a a nd b are real-valued coefficients. The functions 



(fi(x), ipu\(x) and 0i(x) (i — 1, 2) are defined in section 6.1 We then make use of the equations of the WKB 



approximation and the localization assumptions to simplify the expression (7.9). As in the case of fermions, 
this means that the envelopes tpi(x) and (fi(x) are rigidly transported along their respective trajectories and 
so they differ at most by a translation i.e. ipi(x) — (p(x — Xi). 
Again, if we assume the change in the envelope is small 

tp(x) 

for all points x € T>x, the translational difference in the envelopes can be neglected and can be factored out: 

aAf^+bA^x) « <p(x)e ie ^ [a^T,) + b^ 2) (T 2 )e iAe ^) 
where A8 — 02(x) — 6i(x). We then solve the partial differential equation V M 6* = to determine 

6i(x) = k^x* - xt) + Ofa) + 6 (7.10) 

where Oi(xi) — J r k^dx^ are again integration constants. Using that k^ is null and that we are integrating 
along its integral curves, we have J r k^dx* 1 = 0. Thus, the only contribution to the phase difference is 

A6 = k fl (x>{ -x%). (7.11) 

Again note that this phase difference is independent of the position x M € D\ at which the phase difference 
is computed. We also have that the phase difference A9 is independent of the choice of points x\ and x 2 . 
This follows since a change x% — >■ x'/ 1 = x^ + 5x% = x% + €ik^ leaves AO invariant since k^Sx^ = 0. 



As in the fermionic case there is also a phase difference A6>Trans defined by (7.8| related to the transport 



along the trajectories. What is actually measured in a Mach-Zehnder interference experiment is then 

A0 Tot — A9 + A6> Trans . 

Just as in the case for fermions we now neglect the spatial part and we end up with the final qubit quantum 
state at D\ 



I rccomb 



a^ 1) (r 1 ) + ^f 2) (r2)e iAe ^ 



We have now obtained a formalism for describing interference experiments for photons solely in terms of 
two-dimensional quantum states. 

7.2.3 The recipe for adding qubit states 

Above we have established a formalism for quantum interference phenomena for both fermions and photons 
in a Mach-Zehnder interference experiment. This description can be summarized by the following recipe for 
correctly adding the two quantum states: 

1. Transport the quantum states and tp^ to the arbitrary positions x\ and x 2 on the respective paths 



Y\ and T 2 in the recombination region Di using the appropriate transport equation, either (5.11) for 



fermions or (6.11) for photons 



2. determine the integration constants 9\ and 8 2 - For fermions this is determined by equation (7.5). For 
photons this is identically zero. 



3. determine the phase difference A6 using either (7.7) or (7.11) 



4. Finally, the two-dimensional quantum state in region D\ is given by 



/ rccomb 



= a ^ 1} +^ 2) e iAe . 
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7.3 The physical interpretation of phase in terms of wave geometry 

We now provide an intuitive wave-geometric picture for the various terms in the phase difference AO. To do 
this we will focus on the specific case of fermions. Note however that the essential picture is also applicable 



to photons and we will comment on photons when necessary. The phase difference for fermions (7.7) is given 

by 



AO =(fc„ + e A ft) Ax 1 " 1 + f (fc„ + aA M )dx* - f (fc„ + eA^dx" 
Jr 2 Jt! 



= q> k^dx* 1 



< j> A^dx 11 



(7.12) 



where Ax^ = x% — x% and T = T2 + r2_;.i — Ti, where I^-i-i denotes the straight path going from point x^ 
to Xi- The second term in the integral accounts for an Aharonov-Bohm phase. Let us consider the first 
term. The various contributions to this term are 



fc M dx" 



ri 



(7.13) 



The first two terms in the decomposition can each be thought of as representing the accumulation of global 
phase along each trajectory, while the third is related to the displacement of the wavepackets. We now show 
how to interpret these two contributions wave-geometrically. 

7.3.1 The internal phase shift 



The first two terms in ( 7.13 1 are integrals of the wavevector fc M along the paths Ti, i = 1,2. If we parameterize 
the paths with proper time = ^-fc M , the integrals become 



k„dx" 



mc 

dr_ fT 
r n 



(7.14) 



This results in a phase discussed in [211 E2| and motivated from the relativistic path integral. We can also 

a) ti b) t 2 

P. 




Figure 7.2: An illustration of the accumulation of internal phase from a time t\ (a) to a later 
time ti (b) along a trajectory. The internal phase ft nt = Axi nt /R is determined by the difference 
in the offset Axi nt — Ax2 — Axi of a point p p h of constant phase and a point p cnv of constant 
position on the envelope at the two times ti and ti- For timelike packet velocities the phase 
velocity « p h is greater than the group velocity v g = c 2 /v p h < c so the internal phase is seen to 
accumulate along the trajectory. 



understand this term in a simple wave-geometric picture. Consider a point p cnv (T) defined by tp(x) — const 



which is fixed on the rigidly moving envelope (see figure 7.2). With respect to some arbitrary reference frame 
with 4- velocity n' 1 , the velocity at which p cnv moves in this frame is called the group velocity v g (see e.g. 
[50]j. defined by u° = e°u^ = 7„ g = (1 — Vg/c 2 ) -1 / 2 , and corresponds to the particle's velocity. Secondly, 
consider a fixed phase point p p h defined by 6(x) = const. The speed at which this phase point moves is 
given by w p h = c 2 /t> g and is called the phase velocity. Thus, if v g < c the points of constant phase move 
with respect to the wavepacket. It is this difference in velocity that results in the accumulation of the above 
mentioned path integral phase (7.14). To see this, first calculate how much distance Sxi nt is gained by p p h 
relative to p, 



during some time interval dt measured in this reference frame. This is given by 



feint = Pph - Pc 



dt. 
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In order to see how many radians of phase this distance is equivalent to we divide by the reduced wavelength 

X = h/p = h/mv g ; 



SXr 



At 



-7 "^dtrwyVg 



"7 



l dt 



mc 



dr. 



During a finite period of time we have m t = J drmc 2 /% = Ax- m t/X, which is nothing but the path integral 
phase. We can now interpret the path integral phase as how much the constant phase surfaces have shifted 
inside the wavepacket. We call this an internal phase shift mt . When we add two wavepackets it is important 
to keep track of this phase shift as it may lead to destructive or constructive interference. 0- m t calculated for 



each trajectory is simply the integration constants 0i(xi) (e.g. (7.5)) 



Recall that for photons there was no contribution to the phase from the path integral, i.e. the integration 
constants are 0{(xi) = 0. From a wave-geometric picture this is due to the fact that the group and phase 
velocities are equal and therefore mt = . 



7.3.2 The displacement induced phase difference 



Let us now provide a wave-geometric interpretation for the third term in (7.13), fc M Ax M . First, for simplicity 



let Ax^ = Xi — x% be spacelike and orthogonal to some arbitrary unit timclike vector n^, i.e. Ax^ = W v Ax v 



where h?,, = Sti 



projects onto the orthogonal space of n^j^] k^Ax^ then simplifies to: 

X 



\k^Ax^\ = \k fl h' i u Ax v \ = fc^Aa;dis|cos(a) 



where k± — ^—h^ u k^k v and Axd\ s = ^—h^Ax^Ax", and we have used that |cos(a)| = 1 since the 
wavepackets are spatially displaced in the direction of motion, i.e. k^_ oc Ax M . 



Figure 7.3: Illustration of the recombination of two envelopes in the detector region D\. The 
offset of the classical positions xi, X2 of two wavepackets with the same wavelength A produces a 
displacement phase of A#dis = Axdi B /X. 

This contribution to the phase difference, which is present for both fermions and photons, can therefore be 



interpreted as the two wavepackets being spatially displaced, as illustrated in figure 7.3 and as argued in [28] . 
Note that in order for this displacement induced phase 



ference AO^ — k^Ax^ to be experimentally, the 
variance in Aa^dis over runs of an interference experiment must be significantly smaller than the wavelength 
X. Furthermore, if Axdis ~ £ (see Fig |7.1[ ) then the interference effects will be drastically reduced and when 
A^dis > C no interference phenomena will be present. 



7.3.3 Addition of quantum states in the wave-geometric picture 



The recipe for adding two quantum states §7.2.3| can now readily be understood in terms of wave geometry. 
In the detector region Di we have two wavepackets whose envelopes, centred at x± and X2, overlap but 



are slightly offset (as in figure 7.3). Furthermore, each wavepacket has a rapidly oscillating phase that has 
evolved in a path-dependent way along each of the two distinct trajectories Ti and T2 (as in figure 7.2). 
These two effects produce respectively the displacement induced phase difference and the internal phase 
difference. We then add these wavepackets to obtain the total phase difference AO. This is illustrated in 
figure [7T} The total phase difference is again A#Tot = A6>Tr an s. + AO. 



'Ai 11 can always be made spacelike orthogonal by changing the arbitrary end points of the trajectories I?i and Y2- 
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Figure 7.4: From a wave-geometric point of view we can understand the phase difference AO as 
the sum of two contributions: the difference in the internal phase given by A8i nt = Axi U t/R and 
a phase difference A#dis = Axdis/zt originating from the wavepackets being spatially displaced. 



7.4 An example: relativistic neutron interferometry 

As a concrete example for implementing the above recipe for calculating the phase difference we consider the 
gravitational neutron interferometry experiment illustrated in Fig |7.5| known as the Colella-Overhauser- 
Werner (COW) experiment [57]. The setup is geometrically identical to a Mach-Zehnder interferometer: 
The wavepacket is as usual split up into a spatial superposition and the respective wavepackets then travel 
along two distinct paths. The interferometer is oriented such that one path is higher up in the gravitational 
field relative to the other path. Essentially the two components of the spatial superposition have different 
speeds and experience two different gravitational potentials, which leads, in the recombination region, to a 
phase shift. Interference fringes have been observed (see e.g. [27J [251 US] ) when the interferometer is rotated 
in the gravitational field, altering the difference in height of the paths. 





Figure 7.5: A Schematic diagram of a Neutron interferometer used in the COW experiment. 
A neutron incident on the first beam splitter BSi is split into a spatial superposition travelling 
along two distinct paths Ti and rV We find that T2 accumulates a phase shift with respect to 
Ti as it is higher in the gravitational field by Az. 



In this section we are going to derive the phase difference for this experiment using the relativistic 
formalism developed above. The spin of the neutrons is ignored and we treat them as scalar particles with 
no internal discrete degree of freedom. Therefore there is no need to use the transport equation (5.11 ). The 



effects due to the spin could be included by computing (7.8); however the corrections to the overall phase 



difference are minute, as noted in [29) . From this analysis we will arrive at an exact relativistic result which 
contains, in certain limits, both approximate relativistic corrections to the COW experiment [53J [Ml [2j|] as 
well as the non- relativistic result [27J [251 US] ■ 

One might represent the gravitational field for this experiment by the Schwarzschild metric. However, 
since the size of the experimental apparatus is less than a metre and hence small compared to the curvature 
scale, we can mimic gravity by simply going to an accelerated reference frame. This can be achieved by 
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making use of the Rindler coordinates [SDJ 05] in which the fiat spacetime metric takes the forirp^j 



/(1 + ff) 2 



9pv 



\ 






\ 


-1 








-1 





-l) 



Since the Rindler metric is static (i.e. independent of t) we have a Killing vector — (1,0,0,0) 
and hence conserved energy E = p^rf — mc 2 ^goo, where 7 = (500 — fr) -1 ^ 2 with v the speed of the 
neutron as measured in the frame defined by the Killing vector rf . We can now use this conserved energy 
to determine the speed v 2 of the neutron in the upper path given the speed in the lower path v±, i.e. 
Ei = toc 2 7i<7oo(- 2; i) = wc 2 725oo(^2) = E 2 . We take the lower path to be at height Z\ = 0, and so if the 
difference in height is Az we have the height of the top path being z 2 — Az. Therefore we have the relation 
72300(^2) = 7i since goo(zi) = 1. 

The easiest way to calculate the phase difference is by using the formula 



AO 



k^ = 



k^dx" 



k^ 



k^Ax^ 



(7.15) 



where we have used (7.7) and assumed that A^(x) is constant. This formula contains two arbitrary spacetime 
points X\ and x%. Here we take these points to be where the trajectories 1^ and T 2 , respectively, hit the 
second beam splitter BS 2 . Since the spatial positions of these two events are the same in our Rindler 
coordinate system we have Ax M = (At, 0, 0, 0), where At is the difference in arrival time. This is given by 



At = I 



1 1 

Vl v 2 



where £ is the length of the horizontal legs of the paths. The contribution of the phase difference from the 
third term is thus 



mc 2 £ ( 3oo(Az)7 2 g 00 (Az)j 2 



v 2 



\ m,c 2 £ 


( 71 


_ tA 


1 ~ h 




v 2 J 



Let us now turn to the first and second terms in (7.15), representing the internal phase shifts. Since the 



internal phases accumulated along the vertical components of each path are equal the quantity cancels in 
the calculation of the phase difference and so it is unnecessary to calculate them. The internal phase shifts 
of the upper and lower horizontal paths are given by 



g(l) 
7 int 



1 



h h 7ii>i 

The phase difference is then given by 



(2) mc 2 mc 2 t 
int = ~h T2 = "72^2"' 



which simplifies to 



9(1) 

'int 



, . u mc 2 £ f 1 
k u Ax^ = — ' 



1_ 7i _ 7i 

V72^2 71^1 v i v 2 



A9 = 



m^7i 



f'i 



(1 



V2 

_ Azg\2 



(7.16) 



where we can make the replacement v 2 



.9oo 



1 



3oo7i 



Eq.(7.16) is the exact result for the gravitationally induced phase shift in the Rindler metric. This 

i] for the gravitational effect in the COW 
It is instructive to take various limits to 



compares to various results in the literature 



experiment, which turn out to be approximations of (7.16). 
demonstrate these connections. 

5 Note that we could also consider rotating reference frames which would lead to the Sagnac effect |26l 1231 129) . but for 
simplicity we will stick to the Rindler metric. 
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Firstly, if we take the weak field limit, Azg/c 2 <C 1 we obtain 



If we furthermore take Azg/v 2 -C 1, which corresponds to assuming \v\ — t^l/fi <C 1, and consider only 
first order terms, we obtain the phase result in |23j . If we instead take the non-relativistic limit (71 f=a 1) we 
obtain the result of [21] • Expanding this result to two orders of Azg/v 2 gives 

A^?^ + ^V (7.18) 



h \ vi 

which indicates the l g 2 ' correction term derived in 24J. To leading order in Azg/v 2 the non-relativistic limit 
gives the standard theoretical prediction of the phase difference A#cow observed in the COW experiment; 

AW = (7.19) 

There is a reported small discrepancy between measurement and theory [27j . However, the error introduced 
by neglecting corrections in Azg/v 2 is too small to account for this discrepancy [24"] . 



The standard result (7.19) is obtained using a path integral approach (see e.g. [23 [5S]). The path 
integral method allows only for summation over paths which start and end at the same two spacetime 
points. However, the classical trajectories in this problem in fact do not arrive at the second beam splitter 
at the same time. Therefore, the standard expression, although a very good approximation in the specific 
case of the actual experiment under consideration, is not exact even non-relativistically. 



8 Elementary operations and measurement formalism 

In order to develop quantum information theory in curved spacetimes we need to understand how elementary 
operations such as unitary transformations and state updating are represented within the reference frame 
covariant formalism of this paper. This section is dedicated to these issues. In addition we show how 
Hermitian observables are represented, how to calculate their expectation values, and how to construct 
explicitly a quantum observable given the measurement direction of a Stern-Gerlach device, or a polarizer. 



8.1 Fermions 

In this section we develop the notion of unitarity, observables and projectors for fermions. The notion of 
unitarity and observables is not straightforward for two reasons: (1) the inner product is velocity dependent 
and (2) Hilbert spaces associated with distinct points in spacetime must be thought of as separate. The 
formalism that we develop addresses these issues in a reference-frame-covariant way. 



8.1.1 Unitarity and Hermitian operators 



Unitarity is traditionally defined for automorphisms U : H 1— > %, i.e. unitary maps take elements from one 
Hilbert space back to the same Hilbert space. The map U is unitary if it satisfies 



(u<f,\Ufi>) = (m 



(8.1) 



where the inner product is given by (4>\ip) = S A A (f>A'^A- However, for our purposes this definition is too 
restrictive: we are interested in localized qubits transported along some spacetime trajectory T. The Hilbert 
spaces associated with the points along T must be thought of as distinct and the therefore a map induced 



by the transport equation (5.11 1 cannot be thought of as a map from a Hilbert space to itself. Furthermore, 
the inner products for the Hilbert spaces depend on the respective 4-velocity. It is then clear that the 



transformation induced by the transport equation (5.11) is not going to be unitary according to (8.1) as we 
are not dealing with automorphisms. 
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Consider therefore a map U : Hi — > H 2 where Hi and Hi are two Hilbert spaces on the trajectory T. 
That a quantum state belongs to Hi is indicated by a subscript and similarly for Hi- The 'generalized' 
definition of unitarity then becomes 

2 {U<P\Ui>} 2 = i<01V>i (8-2) 

using the inner product l ((t>\tp) l = 4> { A ,I A ^ A i> A , with I A [ A = ufa IA ' A 1 and ,1$ € H u i = 1,2. If we 
adapt the tetrad such that u 1 = (1, 0, 0, 0) along the trajectory we see that the inner product I A A = ui<r IA 
becomes the ordinary inner product 8 A A which is independent of both position and momentum. We would 
therefore expect the inner product between two quantum states along some trajectory V to be conserved. 
To see this let 4>{j) and ip( T ) represent two quantum states that are Fermi-Walker transported, according 



to equation (5.12) along T. Then we have 



Um = v-<w> = uja iA ' A ^-^^ A + „^a'a^,£_^ = (8 _ 3) 



since the Fermi- Walker derivative of u 1 is zero by construction, and the inner product (-|-) is defined using 
Iu(t)- Strictly speaking, the inner product should be labelled with r (i.e. r ("|*)r) m order to indicate that 
we are dealing with different Hilbert spaces. However, for convenience we omit this cumbersome notation. 
Let us now consider a more general evolution dictated by a Schrodinger equation 

D FW ^ A dyj A . fldx» \rIJB, ■ a B 1 

-uj^j + u iaj L IJ A B ip B = lA/tpB (8.4) 



Dt dr \2 dr 

where A J 3 represents some operator on ipg . Requiring the inner product to be preserved under the evolution 
implies that must for all <fi,ip £ H satisfy 

UI (j ia ' a 4>a-A£iI) B - Ul a IA ' A A A B '0 B ,i> A = (8.5) 

or equivalently ((j)\Aip) = (A<j>\ip), which is nothing but the standard definition of a Hermitian operator. 
Since this must hold for all 4>a and ip A we must have I B A A B = I A B A B , . 

In spinor notation we can define A A A = A^l^ B which yields an equivalent definition of Hermiticity 
for the spinorial object A A : 

A A ' A = A A ' A . 

An object A A A satisfying this condition can be written as 

A A ' A = Nja IA ' A (8.6) 

for some real- valued coefficients Ni. We also have A B = I uAA /A A B where I uAA ' is the inverse of I A A 
defined by I uAA iI A B = S A . The corresponding operator A B is then given by 

A B = l uAA ,A A ' B =u I a I AB ,N J a JB ' B =n I N J (^a^ +^>)f 

= -KujNjL 1 -^ B + ujN'S^ (8.7) 

where L IJ ^ are the left-handed s((2,C) generators and the term in 5 B generates changes in global phase 

E3 

It should be noted that the operator A B does not 'look' Hermitian when written out in matrix form. For 
example, A^ ^ A? 1 . Rather, it is only the object A A A = I A B Ag which looks Hermitian in matrix form, 
i.e. A A A — A AA rj\ The reason for this difference can be clearly seen by expressing A B in the rest frame of 
the qubit. In the particle rest frame, Hermitian operators are expressed as A B = Nja° AA ,a IA B , and in this 
case A x 2 = A2 ■ Thus, from an operator A which is Hermitian with respect to S A A we can construct another 



16 Note that one could also have chosen the right-handed representation R ijA b i of the Lorentz group as this would yield the 
same result. 

17 Spinor notation gives the relationship A A ' A = A AA = A. A A between the conjugate and row-column transpose. 
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operator A which is Hermitian with respect to I A A simply by applying a boost, i.e. A B = A A A^ A^ 1B , 
where A is the spin-| representation of the Lorentz boost that takes Sq to u 1 . 

A general inner product preserving evolution can therefore be understood as being composed of two 
pieces. One piece, the Fermi-Walker derivative, dictates how acceleration and the gravitational field affects 
the quantum state and has therefore a purely geometric character. The Fermi-Walker derivative maps 
elements between neighbouring Hilbert spaces %(x,p)(t) an d T~L(x,p)(t+St)- The remaining term A B encodes 
possible non-geometric influences on the quantum state and is an automorphism A : H H. This second 
term is required to be Hermitian with respect to the inner product I A A ^\ 

We have already seen an example of an evolution of the form (|8.4[). In the WKB limit of the minimally 



coupled Dirac equation we arrived at the transport equation (5.11), where the Hermitian operator took on 
the form 

r> B — e R tIJ B _ e h K h L T? T IJB 
d a = -—BjjL A -- — ti I hj t Kh L A 

2m 2m 

where Bjj = hj K h j L F^l is the magnetic field experienced by the particle. To see that the magnetic 
precession term B A B is Hermitian with respect to the inner product, we expand the left side of the Hermiticity 



definition (8.5): 



(4>\6i/)) - = (u K a K B u L IJ - Bui IJ u K a K ) 1>. 

With Bjj real and making use of the identity [a K , L IJ ] = i[r] IJ a K — r] JK a I ] [T], we have 



e 



i^BiP) - (B<t>\V) = ~—u K B u ^[a K ,L IJ ]^ (8.8) 

- -[^^^[ukBj 1 - Bj{uj\ = (8.9) 
2m 

since BjjW 7 — and B/ — 0. The magnetic precession is thus a Hermitian automorphism with respect to 
the inner product I A A . 

8.1.2 Observables and projective measurements 



Observables are represented by Hermitian operators A A , which will take the form indicated in (8.7). The 
covariant expression of the expectation value of the observable A for a spinor %j) A is given by 



ty\A\il>) = ^ AI Nja IA ' A ^ A . (8.10) 



Note that in (8.10) all indices have been contracted, indicating the expectation value is manifestly a Lorentz 
invariant scalar and could in principle represent an empirically accessible quantity. 

In order to complete the measurement formalism we need to discuss how to determine the post-measurement 
quantum state. We do this for the simple case of projection- valued measures, however we can easily extend 
the formalism to generalized measurements. A Hermitian operator has a real eigenvalue spectrum and its 
normalized eigenstates \ip^) are orthogonal, i.e. 

where k, I = ±. The spectral decomposition of an observable A is A A = ^ ± X±P ± A ', where the A± are 
the eigenvalues of A A and the P A represents the corresponding projector onto the eigenstate IV^)- In 
spinor notation, the projectors are given by P^ 3 = I A B 'ip A ,4' A - A pair of projection operators P^ 3 which, 
together with the identity operator, span the space of Hermitian observables on H can also be written as 

P^ = \{5 A B T2mjn J L" A B ) (8.11) 

with njn 1 — 1 and riju 1 = 0. These are the ordinary Bloch sphere projectors but written in a reference 
frame covariant way. One can then suspect that a measurement of spin along some unit direction can be 
represented by such projectors. This is indeed the case as we shall see now in the specific case of Stern- 
Gcrlach measurements. 



18 Hermiticity can alternatively be defined in terms of the partial d/dr or covariant D/Dt derivatives but in doing so we 
would have to modify the definition of a Hermitian operator. 
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8.1.3 The spin operator for a relativistic Stern Gerlach measurement 



To be able to extract empirical predictions from the above formalism we need to determine how N 1 in (8.7) 
corresponds to the relevant parameters defining the experimental setup, e.g. the spatial orientation of a 
Stern-Gerlach magnet. In the literature there exist several proposals for relativistic spin operators and these 
have been studied for various reasons (see e.g. [5H [55J [55J (SH [551 (Ml EES EH]). In this section we are going 
to be concerned exclusively with constructing a spin observable associated with a relativistic Stern-Gerlach 
measurement. Notably, the spin operator that we obtain differs from other proposals and we will elaborate 
on this in a forthcoming paper. 

In order to obtain the correct relativistic spin observable it is necessary to understand in more detail 
the physical aspects of the measurement process. In a Stern-Gerlach spin measurement, a particle is passed 
though an inhomogeneous magnetic field. This causes the wavepacket to separate into two packets of 
orthogonal spin. A subsequent position measurement then records the outcome. To gain further insight, 
let us consider this measurement process in the fermion's rest frame where ej 1 = u^. In such a frame the 
stationary qubit is exposed to a magnetic field B l for a short period of time and it is clear that it is the 
direction (with B 2 = B l B^ij) of the magnetic field that determines what component of the spin we 
are measuring [47] ■ 

If the qubit is moving non-relativistically with respect to the Stern-Gerlach device, the spatial direction 
-g-B 1 of the magnetic field approximately agrees with the orientation of the Stern-Gerlach device m*. How- 
ever, if the qubit is moving relativistically with respect to the apparatus these directions do not necessarily 
coincide, nor is their relationship straightforward. We will establish a relation between these two directions, 
and in doing so we will identify the correct spin observable for a relativistic Stern-Gerlach measurement. In 
particular the relativistic spin operator/observable that we obtain depends on the spatial orientation m 1 of 
the apparatus and the 4- velocities v 1 and u 1 of the apparatus and qubit. 

To proceed we first work out an expression for the electromagnetic field Fjj generated by the Stern- 
Gerlach apparatus. To do that we first introduce the magnetic field 4- vector M 1 = Mm 1 where M is the 
magnitude of the Stern-Gerlach magnetic field. We can now define the electromagnetic tensor as 

Fij = -e LIJK v L M K . (8.12) 



In the rest frame of the apparatus it takes the form 

Fij 





B K 

so there is only a magnetic field, and no electric field, generated by the Stern-Gerlach apparatus in its own 
rest frame. In order to simplify the calculation the gradient e^V^M of the magnetic field strength is assumed 
to point in the same direction as the magnetic field itself, i.e. e^V^M oc mjQ 

We can now calculate the magnetic field 4-vector B 1 corresponding to the magnetic field as measured in 



the rest frame of the qubit. It is given by B 1 = \e LIJK ulFjk ■ Using (8.12), we obtain 



B 1 = - l -e LIJK u L e MJKN v M M N = M 1 {v ■ u) - v 1 (M ■ u) 

with e LIJK e MJKN — —2(5^5^ — Sn^m) and e 0123 = 1 H5J p87]. For spin measurements, the 4-vector n 1 is 
now the normalized qubit rest-frame magnetic field, i.e. 

it \ — B ' 

n {m,u,v) = — 

where B = \J — B I B ,J r\ 11 ^ and it is easy to check that n • u = 0. This expression n 1 becomes singular only 
for unphysical or trivial situations characterized by u 1 being null, or M = 0. 

Thus, given the spatial orientation of the Stern-Gerlach apparatus and the 4-velocities of the apparatus 



and qubit we obtain, using (8.7), a relativistic spin operator given by 



S A B = -2iu I n J (m,u,v)L IJ A B . 



The expectation values are calculated using (8.10) and the corresponding projectors are given by (8.11). 
We now have a fully relativistic and reference frame invariant measurement formalism for a Stern-Gerlach 
measurement. 



19 Although this is not strictly possible as the magnetic field must satisfy Vim 1 = one can always choose a field which 
approximately has eJ'VjjM oc mj locally 1621 . 
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8.2 Photons 

In this section we develop the notion of unitarity, observables and projectors for photons. The definition 
of unitarity is in some sense simpler than for fermions as the inner product is not velocity dependent; the 
difficulty is only in handling the gauge degrees of freedom. 

As we have already stated, the polarization state of a photon can be represented by a spatial complex 
4- vector ip 1 orthogonal to the null wavevector . Furthermore, the corresponding quantum state of a pho- 
tonic qubit with null velocity u 1 was identified as being a member of an equivalence class of polarization 



vectors ip 1 ~ ip 1 + vu 1 all orthogonal to u 1 (\ 6.4 1. With the orthogonality condition and the gauge degree of 
freedom this space therefore reduced to a two-dimensional Hilbert space on which unitary and Hermitian op- 
erators act. We now develop the notion of unitarity, observables and projectors within this four-dimensional 
formalism. 



8.2.1 Unitarity and Hermitian operators 



Unitarity and Hermiticity are more straightforward with polarization vectors than spinors because the defi- 
nition of unitarity (U<p\Uip) = (4>\ip) is in terms of a standard inner product rjjj^ip" 1 where rjjj is constant. 
The requirement for unitarity again translates into requiring that the inner product between two polarization 
vectors is conserved along trajectories, i.e. 



dA DX DX 

Consider now a Schrodinger evolution of the form 



Dip 1 
~DX 



0. 



Dip 1 

Jdx 



= Pu 1 + iA J jip J 



(8.13) 



(8.14) 



which is more general than the transport equation (6.111. Substituting into (8.13), we get 

A/^jiP 1 



rA 1 ^ = 



where we have used that ip^j — u 1 (pi = 0. Requiring that this hold for all (p 1 and ip 1 we obtain the standard 
definition of Hermiticity; A/ — A 1 j. 

In addition to the above we also require that the gauge condition uj-ijj 1 = be preserved. This implies 
that 



d(u/i//) Dimip 1 ) Dip 1 



dA 



DA 



DA 



ui = A I J ip J u I = 0. 



This condition ensures that Hermitian operators A j map polarization vectors into polarization vectors. 
Again this should hold for all ip J , so we have the condition A 1 jUj oc uj. The following two conditions suffice 
for characterizing a general Hermitian operator; 



Aj — A j, 



A 1 



u oc u 1 . 



.15a) 
.15b) 



In order to determine the form of valid operators it is convenient to express the matrix A I J in terms of 
w I ,f(,fj} which spans the full tangent space. Recall from Section 



6.4.2 



that a 



a set of basis vectors {u 

diad frame f J A = (// , f£) defines a spacelike two-dimensional subspace orthogonal to two null vectors u 1 , w 1 
with — 1. f^fj = hj = $j~ ufwj — w I uj is the metric on the spacelike subspace |51j . One can then 

define a sixteen-dimensional complex vector space spanned by the outer products of {u 1 ,1V 1 , // , f£} with 
the dual vectors {uj,wj, fu, /aj}. Components of an arbitrary matrix A I J G (B ® B) of this space can 
then be expanded in terms of these sixteen elements. 

A valid map A I J on polarization vectors must satisfy equations 



.15). In terms of the sixteen basis 



elements, no terms in w 1 or wj can exist, since w I ki 7^ 0. Any remaining terms that involve u 1 or uj 
are pure gauge and do not change the polarization vector. A hermitian operator is therefore, up to gauge, 
represented as 

A 1 ., = aflfj + pflft + Bflfj + bflfi 
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where the real numbers a and b and the complex number j3 constitute the four remaining real degrees of 



freedom. The two conditions in (8.15) thus reduce a 4 x 4 hermitian matrix to effectively a 2 x 2 hcrmitian 



operator that acts on the transverse spacelike (polarization) degrees of freedom. Such an operator can be 
written in terms of the Pauli matrix basis as A A B — C a a a A B , where the Pauli matrices act on the two 
dimensional Jones vector. C a consisting of four coefficients a = 0, 1, 2,3 does not transform as a 4-vector 
and therefore does not have any spatial significance. For the operator where a = 1, 2 or 3 the eigenbasis 
corresponds to respectively the diagonal linear polarization basis, the circular polarization basis, and the 
horizontal-vertical polarization basis of the Jones vector. The relation between the four-dimensional and 
two-dimensional hermitian operators is A I J = fA_A A B ff. 

8.2.2 Observables and projective measurements 

The construction of observables and projectors is identical to that of fermions. Observables are represented 
by Hermitian operators A 1 j. Let Pl^ represent the eigenvectors of A 1 j. The P^ form an orthonormal basis 
with Plj^P^j — Ski- The probability p^ of getting outcome X k is given in tetrad notation by 

Pk = \P( k) M 2 - (8-16) 

The corresponding projector for an eigenvector PL^ is Pj = PhsPy^j and the post-measurement state, up 

to gauge, is given by if) 1 —> ip' 1 = P I J ip J . 

In this case we have a clean interpretation of projectors Pj. as polarizer vectors P 1 : complex, spacelike 
normalized vectors orthogonal to photon velocity, P I uj = 0. Polarizer vectors correspond to the physical 
direction and parameters of an optical polarizer: A linear polarizer direction is of the form e lfl P 1 with P 1 
real, and a circular polarizer is a complex vector P 1 = ^(P± + i-P/) with PiP^Vu — and PfPfrju = 
P2P2VU = 1- The probability of transmission of a polarization vector through a polarizer is simply the 



modulus square of the overlap of the polarization state with the polarizer vector (8.16). Such an overlap 
clearly does not depend on the tetrad frame used, and indeed all tetrad indices are contracted in (8.16). 
The probability p is then manifestly a Lorentz scalar. It is easy to verify that the formalism is invariant 
under gauge transformations ip 1 — > ip 1 + vu 1 and P 1 — > P 1 + ku 1 . Thus, the probability p is both gauge 
invariant and Lorentz invariant as should be the case. With this completed measurement formalism it is 
then not necessary to work with the Jones vector or Wigner rotations, both of which involve working with 
the cumbersome adaption procedure. 



9 Quantum entanglement 

Until now we have been concerned with the question of how the quantum state of some specific physical 
realization of a single qubit is altered by moving along some well-defined path in spacetime. We shall now 
show how this formalism can easily be extended to describe entanglement of multiple qubits. 

9.1 Bipartite states 

We have seen that it is necessary to associate a separate Hilbert space with each pair of position and 
momentum (x^^p^). A single qubit moving along a specific path x^(A) in spacetime will therefore have its 
state encoded in a sequence of distinct Hilbert spaces associated with the spacetime points along the path 
The formalism that we have so far developed determines how to assign a quantum state to each distinct 
Hilbert space along the path along which we move the qubit. The one-parameter family of quantum states 
|^(A)) is parameterized by some parameter A of the path x(A), and the sequence of Hilbert spaces associated 
with the path is H(x,p)(X)- The quantum state \ip(\)) belongs to the specific Hilbert space H( x ,p)(\)- 

Let us now consider how to generalize the formalism of this paper to the quantum state of two, possibly 
entangled, qubits in curved spacetime. Instead of one worldlinc we will now have two worldlines, Xi(Ai) and 
£2(^2), and consequently instead of one parameter A we now have two, Ai and A2. Corresponding to each 

20 Recall that while we can uniquely determine the 4-momentum = m^^P from the trajectory x(t) in the case of massive 
fermions, the same is not true for photons. Due to the arbitrariness of the parametrization of the null trajectory x(X) we can 
only determine the null momentum up to a proportionality factor, i.e. p M <x . 
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value of Ai and A2 we have two spacetime points and two Hilbert spaces, %(xi,pO(Ai) an d H(x 2 ,p 2 ){\ 2 )- ^ ^ s 
therefore clear that the quantum state describing the two qubits is mathematically described by a quantum 
state \4>{\i, A2)) £ T~L{x 1 ,p 1 ){\ 1 ) ® H{x 2 ,p 2 ){\ 2 ) which belongs to the tensor product Hilbert space. 

In order to calculate statistics (e.g. correlation functions) we need to provide an inner product for the 
tensor product Hilbert space %( Xl , Pl )(\ 1 ) <£> '^L(x 2 ,p 2 )(x 2 )- The natural choice is the inner product induced by 
the inner products for the individual Hilbert spaces. 

In the case of fermions the natural choice for the parameterization A is the proper time r. The Hilbert 
spaces 'H(x 1 .p 1 )(T 1 ) an d H(x 2 , P2 )(t 2 ) nave tne mner products given by I A A = u\a IA A and I A A — u\a IA A 
where u\ and uj are the respective 4-velocities. In our index notation a bipartite quantum state can be 
represented by an object with two spinor indices iPa 1 a 2 ( x i{ t i)t x 2{i~2))- Note however that the indices A\ 
and Ai relate to two distinct spinor spaces associated with two distinct points x\(t\) and xiij^) and therefore 
cannot be contracted. The inner product between two bipartite quantum states l^) and |</>) becomes 

(M4>)puP2 = u}a IA '^u 2 ja JB * B *Tp A , lB ^ AlB2 . 

where pi, P2 are the momenta of the two qubits. In the case of photons the quantum state \tp) can be 
represented by a polarization 4- vector ?/> 7 (A). A bipartite state is then given by a two-index object 
?/> 7l/2 (a;i(Ai), £2^2)), where I\ and I2 belong to two different tangent spaces and thus cannot be contracted. 
The requirement that the polarization vector be orthogonal to the null wavevector generalizes to u^ifj 1112 = 
= uj 2 ip I±l2 . The inner product between two bipartite quantum states and \4>) becomes 

(M4>)vuP2='nhJ 1 m 2 j,3 Ill2 4> JlJ2 - 



We could also consider bipartite states \<p) where one component is an electron and the other is a photon. 
Mathematically this would be represented as <fy (xi(\),X2(t)) and the inner product can be constructed 
similarly. 

Let us now turn to the evolution of bipartite quantum states. The physically available interactions of 
the qubits are given by local operations. Mathematically this means that the most general evolution of the 
state vector is given by two separate Schrodinger equations: 

i^H^(Ai,A 2 )) - li(A 1 )®I^(A 1 ,A 2 )) (9.1) 
i^|V(Ai,A 2 )) - I®i 2 (AaMAi,A2)) (9- 2 ) 



where Ai(Xi) and A 2 (X2) are possible local Hermitian operators (as defined in |8j) acting on the Hilbert 
spaces Ti(x 1 ,p 2 )(Xi) an d ^-(x 2 .p 2 )(x 2 )- D T /DX denotes the transport law, i.e. the Fermi-Walker transport for 
fermions or the parallel transport for photons. 

This mathematical description of the evolution of the quantum state is not standard since we have 
two Schrodinger equations rather than one. However, if we introduce an arbitrary foliation t{x) these two 
equations can be combined into one Schrodinger equation. First we express the parameters as functions of 
the foliation Ai = X\(t) and A2 = X2{t). This allows us to write the quantum state as only depending on one 
time parameter: \ip(t)) = \tp(Xi(t), A 2 (t))). The evolution of the quantum state now takes a more familiar 
form 

■ D T , ,. N . . /dAi D T dA 2 D T \ . ,„ , . , . ... 

^LA 1 <g> I + ^1 ® A^j \yj(X 1 (t),X 2 (t))} 

where A = ^j-Ai®I+ ^-I(S)j4 2 is the total Hamiltonian acting on the full state. For this single Schrodinger 
equation to hold for all paths a^(Ai) and £E 2 (A 2 ), and all choices of foliation t(x), and so for all values of ^± 



and it is necessary that both equations (9.1| and (9.2) hold. Thus, the two mathematical descriptions 



of the evolution of the quantum state are equivalent when only local operations enter in the evolution. 

If the Hamiltonian is not a local one, i.e. not of the form A = aA\ ® I + bl (g> A 2 , then it is not possible to 
cast it into the previous form with two independent evolution equations and it is also necessary to introduce 
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a preferred foliation. However, if all interactions are local the introduction of an arbitrary foliation is not 
necessary. 

The generalization to multipartite states is straightforward. Furthermore, if we are dealing with identical 
particles the wavefunction should be symmetrized or antisymmetrized with respect to the particle label, 
depending on whether we are dealing with bosons or fermions. This will correctly reproduce the Pauli 
exclusion phenomenon and the Hong-Ou-Mandel bunching phenomenon |63j . 



9.2 State updating and the absence of simultaneity 

Let us now discuss the issue of state updating for entangled states. Let Ti and T 2 be two spacetimc 
trajectories along which two qubits are being transported. Furthermore let xi(Ai) and ^2(^2) each represent 
a distinct point on the corresponding trajectory. Consider now that the two qubits are in the entangled state 

l^(Ai,A 2 )) = (l+.AO |-,A 2 ) - I- Ai> |+,A 2 )) . 

If a measurement on qubit 1 is carried out at the spacetime point xi(Xi) with outcome '+', the bipartite 
state has to be updated as follows: 

|V(Ai,A 2 )) |+,Ai)|-,A 2 ). 

Note that the value of A 2 is left completely arbitrary after the state update. It should therefore be clear 
that even though the state update is associated with a distinct spacetime point xi(Ai) on the trajectory Fi 
no such point can be identified for T 2 . In other words, as long as the local unitary evolution for particle 2 
is well-defined, the state updating can be thought of as occurring at any point £2^2) along T2- The point 
£2^2) could be in the past, elsewhere, or in the future of Xi(Ai). The reason for this freedom in state 
updating is that a projection operator on particle 1 commutes with any local unitary operator acting on 
particle 2. 



9.3 An example: quantum teleportation 

As an example of entanglement and state updating let us look at quantum teleportation in a curved spacetime. 
Although the mathematics is virtually the same as for the non-relativistic treatment, the interpretation is 
more delicate. In particular, in a curved spacetime the claim that the input state is in some sense the "same" 
as the output state seems to lack a well-defined mathematical meaning. However, as we are going to see, 
in order to carry out the standard teleportation protocol the parties involved must first establish a shared 
basis in which the entangled state takes on a definite and known form. For example, Alice and Bob could 
choose the singlet state. This will be called the 'canonical' form of the entangled state. Once this shared 
basis has been established, Alice and Bob have a well-defined convention for comparing quantum states 
associated with these different points in spacetime. The problems associated with quantum teleportation in 
curved spacetime are therefore similar to the problems associated with teleportation in flat spacetime when 
the maximally entangled state is unknown [BH ESI EE1 EZ] • 

Consider then three qubits moving along three distinct trajectories Ti, T2, and T3. For concreteness 
assume that the qubits are physically realized as the spins of massive fermions. The tripartite state is then 
given by 

|T; Ai, A 2 , A 3 ) G %(xi,pi)(Ai) ® K(x 2 , P2 )(\ 2 ) ® ^(x 3 , P3 )(x 3 ) 

or, written in our index notation, T a 1 a 2 a 3 (^i, A 2 , A3). 

In order to proceed we define a basis for each one of the three Hilbert spaces 'H( Xl , Pl )(x 1 ), T~L{x 2 ,p 2 ){\ 2 )^ 
and %(x 3 ,p 3 )(A 3 ) an d for all points along the trajectories Ti, T2, and T3. The three pairs of basis vectors are 
assumed to be orthonormal, and to evolve according to the local unitary evolution (e.g. by pure gravitational 



evolution given by the Fermi- Walker transport ( 5.12 1). Therefore, once we have fixed the basis for particle 



i at one point iCj(Aj) on the trajectory T,, the basis is uniquely fixed everywhere else along the trajectory, 
at least where the local unitary evolution is well-defined. It follows that a state can be expressed as a linear 
combination of these basis states with the components independent of Aj, i = 1, 2, 3. 
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We denote these three pairs of (one-parameter families of) orthonormal states as </>^(Ai) and (A$), 
with i = 1, 2, 3. The orthonormality conditions are explicitly given by 

(0«;A^«;A 4 ) = «» ^^(A^g (A,) = 1 
(V^Aih^jA,) = t»i* ) ^ A< ^2j ) (W2(A i ) = l 

with i = 1,2, 3. 

Consider now the specific tripartite three-parameter family of states 



r^WAi.Aa.Aa) = 7s (Wi^Ai) + M^Ai)) (^ 2 ) (A 2 )0g(A 3 ) + ^(A^^^As)) (9.3) 



where a and (3 are independent of Ai (since evolution is entirely in the basis vectors) and \a\ + \(3\ = 1. 
The maximally entangled state (involving particles 2 and 3) only has coefficients equal to or 1 which are 
trivially independent of the parameters A2 and A3. We will call the maximally entangled state in equation 



(9.3) the canonical form. Other choices of this canonical form are possible but the teleportation protocol 



used below (see (9.4)) will then change accordingly. 

As in the flat spacetime description we now proceed to rewrite the state in the Bell basis for the Hilbert 
space H( X1)P1 )(A 1 ) ® ft(x a ,p a )(A 2 ) : 



*1a 2 (Ai,A 2 ) ee ^(^ 1 1 ) (A 1 )^(A 2 )±< ) (A 1 )^2(A 2 ) 
*1a 2 (Ai,A 2 ) ee ^(^(Ax^gCA^i^CWSCA,)). 
In these new bases the state |T) reads 

Ta 1 A 2 A 3 (Ai, A 2 , A3) = ^($J lA2 («0A 3 +WA3) + *A!A 2 W'Aa ~ Pi>A 3 ) 

+ ^Xa 2 {0<Pa* + a^ A3 ) + {-P4>a, + a^ As )). 

Alice now performs a Bell basis measurement on the particles 1 and 2. It is important to note that the 
specific physical measurement operation that Alice needs to carry out depends on the bases (<Pa^^a^) anc ^ 

(2) (2) (2) (2) 

(4> A ^,ip A ^). The basis {4>\ 2 , wa^) is determined by the maximally entangled state. If Alice does not know 
the maximally entangled state she will not be able to do the correct Bell basis measurement. 

The outcome of the Bell basis measurement ($ + , V I' + , or is then communicated to Bob's side: 
Bob performs the local unitary operation U (which is assumed to act on the state but not on the basis) 
given by 




U = 7 (9-4) 



where a x ,a y , a z take the usual form when expressed in the local orthonormal basis {4>aI > ^aI)- Thus, as was 
the case for Alice, in order for Bob to know which specific physical operation to carry out, the basis has to 
be specified. Since the basis is determined by the maximally entangled state, Bob must know the entangled 
state in order to carry out his operations. Note that if the basis in which Bob applies the operation is 
incorrect, the state that Bob obtains at the end differs depending on the outcome of the Bell measurement. 
In the case where Bob implements the operation in the correct basis, Bob's state is given by a4>A 3 + fi^A?,- 
We should now ask whether we can sensibly view this protocol as a 'teleportation' of a quantum state 
from Alice to Bob in the sense that Bob received the same state as Alice sent. This hinges on there being 
a meaningful way of comparing quantum states associated with distinct spacetime points. However, as we 
have already stressed, if spacetime is curved, sameness of quantum states cannot be established uniquely 
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by parallel transporting one qubit to the other as this would depend on the specific path along which we 
transport the qubit. 

On the other hand, the maximally entangled state, by determining the bases for it to take the canonical 
form, defines a shared spinor basis for Alice and Bob. If we change the basis, the maximally entangled 
state would of course change accordingly and would no longer take on the canonical form. Given a shared 
basis we have a well-defined way of comparing quantum states and in particular a well-defined way to claim 
that they are the same or not. Thus, when we have a maximally entangled state, there is a natural way 
of comparing quantum states associated with distinct spacetime points. It is in using this convention for 
comparing quantum states that we can claim that Bob did indeed receive the same quantum state, and 
therefore we can say that the state was in this sense teleportedp"| 

In the case of fermions it is also easy to see that the maximally entangled state will also establish a 
shared reference frame, i.e. a shared tetrad. This comes about because from the left-handed spinor ipA by 
means of which the quantum state is expressed we can construct the null Block ^-vector b 1 — a IA A, tpA''4 , A- 
A maximally entangled state can therefore be loosely understood geometrically as a kind of 'non-local 
connection'. 



10 Conclusion, discussion, and outlook 

Recently there has been increased interest in exploring relativistic quantum information theory in the context 
of phenomena from quantum field theory such as the Unruh effect and particle number ambiguity [6l [9] . In 
contrast, this paper explored relativistic quantum information in the regime where such effects are negligible 
and restricted attention to localized qubits for which the particle number ambiguity is circumvented. A 
localized qubit is understood in this paper to be any object that can effectively be described by a position 
and momentum (x,p) and some two-component quantum state \ip). We obtained a description of localized 
qubits in curved spacetimes, with the qubits physically realized as the spin of a massive fermion, and the 
polarization of a photon. 

The original motivation for this research was to develop a formalism for answering a simple experimental 
question: if we move a spatially localized qubit, initially in a state l^i) at spacetime point x\, along some 
classical spacetime path T to another point x%, what will the final quantum state be? Rather than 
working directly with Wigner representations our starting point in answering this question was the one- 
particle excitations of the quantum fields that describe these physical systems. The one-particle excitations 
in curved spacetime satisfy respectively the Dirac equation minimally coupled to the electromagnetic field, 
and Maxwell's equations in vacuum. From these fields we were able to isolate a two-component quantum 
state and a corresponding Hilbert space. 

In the case of fermions, the equation governing the transport of the spin of a fermion consisted of a spin- 
5 version of the Fermi-Walker derivative and a magnetic precession term, expressed in a non-orthonormal 
Hilbert space basis. This result was expected since an electron can be regarded as a spin-| gyroscope, and 
the precession of a classical gyroscope along accelerated trajectories obeys the Fermi-Walker equation. By 
introducing an orthonormal Hilbert space basis, which physically corresponds to representing the spinor 
in the particle's rest frame, we reproduced the transport equation obtained in [101 111] which made use of 
Wigner representations. 

We showed by applying the WKB approximation to vacuum Maxwell equations that the polarization 
vector of a photon is parallel transported along geodesies and that this corresponds to a Wigner rotation. 
Furthermore, this rotation is proportional to the spin-1 connection term uinii when we consider a reference 
frame where the photon 3-velocity is along the z-axis. In this way the effect of spacetime geometry on the 
quantum state was easily identified. 

We worked with faithful finite-dimensional but non-unitary representations of the Lorentz group, specif- 
ically a two-component left-handed spinor ipA, and polarization 4- vector tjj 1 . Nevertheless, by identifying 
a suitable inner product we obtained a unitary quantum formalism. The advantage of working with non- 
unitary representations is that the objects which encode the quantum state transform covariantly under 
actions of the Lorentz group. As a result the transport equations are manifestly Lorentz covariant, in ad- 
dition to taking on a simple form. The connection to the Wigner formalism was obtained by choosing a 



We note that all examples of experimentally produced entangled qubit pairs are produced in localized spatial regions and 
distributed to the parties. The components of the entangled state will then undergo local unitary evolution along each trajectory. 
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reference frame adapted to the particle's 4-velocity, which is a reason why the Wigner formalism is not 
manifestly Lorentz covariant. 

In order to make empirical predictions we need a way to extract probabilities for outcomes from the 
formalism. Such a measurement formalism was developed for both fermions and photons. The predicted 
probabilities of outcomes of experiments were shown to be manifestly Lorentz invariant and thus reference 
frame invariant, resulting in a relativistically invariant measurement formalism. We also derived the specific 
Hermitian operator corresponding to a Stern-Gerlach measurement, by providing a physical model of this 
measurement process. In this way a unique spin operator can be identified given the spatial orientation and 
velocity of the Stern-Gerlach apparatus and the velocity of the particle. Notably this operator does not 
agree with previous competing proposals [59, 60, 61 . 

A second advantage of working in terms of the Dirac and Maxwell fields instead of the Wigner representa- 
tions is that global phases and quantum interference come out automatically from the WKB approximation. 
By considering spacetime Mach-Zehnder interference experiment we arrived at a general relativistic formula 
for calculating the gravitationally induced phase difference. In the specific case of gravitational neutron in- 
terferometry we reproduced the existing formulae for the gravitationally induced phase difference as various 
limits of our formula. Our overall approach, however, provides a general, unified, and straightforward way 
of calculating phases and interference for any situation. 

Finally we generalized this formalism to the treatment of multipartite states, entanglement, and telepor- 
tation, thereby extending the formalism to include all the basic elements of quantum information theory. 

The Lorentz group played a primary role in the construction of qubits in curved spacetime. This role 
can be understood in terms of how gravity acts on physical objects. When an object is moved along 
some path x ,J '(X) in spacetime it passes through a sequence of tangent spaces. These are connected by 
infinitesimal Lorentz transformations that are determined from the trajectory and the gravitational field 
(i.e. the connection 1-form u>J 7 ). That is, apart from a possible global phase, this sequence of infinitesimal 
Lorentz transformations determines how an object is affected by the gravitational field. In particular, if the 
object has internal degrees of freedom that transform under the Lorentz group, we can determine the effect 
of gravity on the state of these internal degrees of freedom. For example, this is the explanation for the 
presence of spin connection terms in the fermion Fermi- Walker transport. It is also an explanation for why 
the photon Wigner rotation is simply a rotation of the linear polarization and so respects the helicity of the 
photon: no Lorentz boost can change frames sufficiently to change the helicity of a photon. 

In this paper we focused on just two physical realizations which constituted non-trivial representations of 
the Lorentz group. We can nevertheless contemplate other realizations such as composite two-level systems. 
In order to understand how gravity acts on the qubit state the same general approach applies: One needs 
to provide a mathematical model of the physical system. Once a model is established one can in principle 
determine how (if at all) the quantum state transforms under a Lorentz transformation. In addition to this 
there are other possible gravitational influences on the quantum state such as gravitationally induced phases. 

As a concrete example of a physical realization that would behave very differently to the elementary 
realizations treated in this paper, consider a two-level system where the two levels are energy eigenstates \E\) 
and | E 2 ). From ordinary non-relativistic quantum mechanics we know that the total state \ip) = a \Ei)+b \E 2 ) 
will undergo the evolution 

|V(*)> = ae iElt/h \E X ) + be lE ^l h \E 2 ) . 

If the composite object is much smaller than the curvature scale and the acceleration is sufficiently gentle 
to not destroy it we can obtain a fully general relativistic generalization by simply replacing the Newtonian 
time t with the proper time r: 

|^(r)) = ae iElT/h \Et) + be iE2T/h \E 2 ) . 

Thus, because of the energy difference we develop a relative phase between the two energy levels which is 
proportional to the proper time of the trajectory. In principle we can make use of such a two-level system 
to measure the proper time of a spacetime trajectory. Since proper time is path dependent we see that the 
transport of the quantum state is also path dependent, and we can also contemplate possible interference 
experiments. 

The formalism presented in this paper provides a basis for quantum information theory of localized qubits 
in curved spacetime. One theoretical application of this is to extend the applicability of clock synchronization 
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[551 IM1 170] and reference frame sharing [65, 66 to curved spacetime, and provide an interesting physical 
scenario for the study of symmetries in quantum mechanics [TTJ [72] . 

This formalism also has potential measurement applications. For example, in the last ten or fifteen 
years there has been interest in precision measurement of the effects of general relativity. Most recently 
there has been the experimental confirmation of the predicted frame-dragging effect by Gravity Probe B 
[73 Q On the other hand, in the same period there has been an increased interest in quantum precision 
measurements using techniques from quantum information theory [751 1761 1771 1781 1791 180) . The formalism 
of this paper provides a bridge between these developments, providing a solid foundation for considering 
the effects of gravity on quantum states, and for considering the design of precision measurements of these 
effects. Importantly, by the use of entangled states one can in principle significantly increase the precision 
of such a measurement. For localized qubits in curved spacetimes as defined in this paper, the effect of 
gravity enters as classical parameters in the unitary evolution of the quantum state. Therefore, one should 
be able to use these same quantum information theory techniques to increase the precision in measurements 
of the gravitational field. It is plausible that such an amalgamation of the transformation of the discrete 
degrees of freedom of a quantum state and the phase accumulation, by increasing the degrees of freedom 
to be measured, will increase the sensitivity with which possible future sophisticated precision quantum 
measurements can measure gravitational effects. 

For example, spacetime torsion is generally believed, even if non-zero, to be too small to measure with 
present day empirical methods [45 . One problem is that torsion, as it is conventionally introduced in 
Einstein-Cartan theory, does not have any propagating degrees of freedom. Thus, the torsion in a spacetime 
region is non-zero if and only if the spin density is non-zero there. Experiments to measure torsion thus require 
objects to pass though a material with non-zero spin density to accumulate an effect, while accounting for 
standard interactions. Needless to say, measuring torsion is then very difficult. However, electrons decouple 
from matter in the high energy WKB limit and effectively only feel the gravitational field including torsion. 
Thus, the spin of high energy fermions might carry information about the spacetime torsion. 

Finally, on the more speculative side it would perhaps also be interesting to model closed timelike curves 
[81] l82l [83] within this formalism. It is possible that the relativistic approach to state evolution and bipartite 
states provides rules or constraints for the manipulation of quantum information in closed timelike curves. 

To summarize: in this paper we have provided a complete account of the transport and measurement 
of localized qubits, realized as elementary fermions or photons, in curved spacetime. The manifest Lorentz 
covariance of the formalism allows for a relativistic treatment of qubits, with a perhaps more straightforward 
interpretation than approaches based on Wigner representations. The treatment of multipartite states, 
entanglement and interferometry provides a basis for quantum information theory of localized qubits in 
curved spacetime. 
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Appendix 

A Spinors and 5L(2,C) 

In our analysis of qubits in curved spacetime it will be necessary to introduce some notation for describing 
spinors. A spinor is a two-component complex vector (f>A, where A — 1,2 labels the spinor components, living 
in a two-dimensional complex vector space W. We are going to be using spinors as objects that transform 
under SL(2,C), which forms a double cover of SO + (l,3). Hence, W carries a spin-i representation of the 
Lorentz group. The treatment of spinors in this section begins abstractly, and ends with details specific to 
Dirac spinors. 



See also 1741 1131 for a recent theoretical analysis of polarization rotation due to gravity. 
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A.l Complex vector spaces 

Mathematically, spinors are vectors in a complex two-dimensional vector space W . We denote elements of 
W by 4>a- Just as in the case of tangent vectors in differential geometry, we can consider the space W* of 
linear functions ip : W n- C, i.e. ip(a4>i + (3 $2) = &4>(4>i) + P^ifa)- Objects belonging to W* , which is 
called the dual space of W, is written with the index as a superscript, i.e. i> A G W* . 

Since our vector space is a complex vector space it is also possible to consider the space W of all 
antilinear maps x '■ W 1— > C, i.e. all maps x such that x( a 0i + P<j>2) — oOii^i) + Pxifa)- A member of that 
space, called the conjugate dual space of W, is written as x G W . The prime on the index distinguishes 
these vectors from the dual vectors. 

Finally we can consider the space W dual to W , which is identified as the conjugate space of W. Members 
of this space are denoted as £4/ . 

In summary, because we are dealing with a complex vector space in quantum mechanics rather than a 
real one as in ordinary differential geometry we have four rather than two spaces: 

• the space W itself: <f)A G W\ 

• the space W* dual to W: ip A eW*; 

• the space W* conjugate dual to W: x A S W*; 

• the space W dual to W* : £4/ G W. 



A. 1.1 Spinor index manipulation 

There are several rules regarding the various spinor manipulations that are required when considering spinors 
in spacetime. Specifically, we would like to mathematically represent the operations of complex conjugation, 
summing indices, and raising and lowering indices. The operation of raising and lowering indices will require 
additional structure which we will address later. 

Firstly the operation of complex conjugation: In spinor notion the operation of complex conjugation will 
turn a vector in W into a vector in W. The complex conjugation of <pA is represented as 

<t>A = ~4>A> ■ 

We will also need to know how to contract two indices. We can only contract when one index appears 
as a superscript and the other as a subscript, and only when the indices are either both primed or both 
unprimed, i.e. 4>a?P A and £,A'X A are allowed contractions. Contraction of a primed index with an unprimed 
one, e.g. 4>aX A j is n °t allowed. 

The reader familiar with two-component spinors [841 [Tl 155] will recognize the index notation (with primed 
or unprimed indices) presented is commonly used in treatments of spinors. It should be noted however that 
this structure has little to do with the Lorentz group or its universal covering group SX(2,C). Rather, this 
structure is there as soon as we are dealing with complex vector spaces and is unrelated to what kind of 
symmetry group we are considering. We will now consider the symmetry given by the Lorentz group. 



A. 2 SL(2, C) and the spin-^ Lorentz group 

The Lie group SL(2, C) is defined to consist of 2 x 2 complex-valued matrices Cj^ with unit determinant 
which mathematically translates into 

1 ab r c r d 1 

-eerie L> A L B - 1 

where e AB is the antisymmetric Levi-Civita symbol defined by e 12 = 1 and e AB = —e BA and similarly for 
sab- It follows immediately from the definition of SL(2,C) that the Levi-Civita symbol is invariant under 
actions of this group. If we use the Levi-Civita symbols to raise and lower indices it is important due to 
their antisymmetry to stick to a certain convention, more precisely: whether we raise with the first or second 
index. See e.g. [T] or [S3] for competing conventions. 

The generators G IJ in the corresponding Lie algebra s((2, C) is defined by (matrix indices suppressed) 

[G /J , G KL ] = i (r, JK G IL - V IK G JL ~ V JL G IK + V IL G JK ) 
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and coincides with the Lorentz so(l,3) algebra. In fact, «SX(2,C) is the double cover of SO + (l,3) and 
is therefore a spin-| representation of the Lorentz group. Note also that the indices /, J,K,L = 0, 1,2,3 
labelling the generators of the group are in fact tetrad indices. The Dirac 4x4 representation of this algebra 
is given by 

where the r y I are the usual 4x4 Dirac 7-matrices. This representation is reducible, which can easily be seen 
if we make use of the Weyl representation of the Dirac matrices 

r_( g\ a , 
7 - I -I A A 



in which the generators become 



4 L' " J ^ (R IJ ^ A ' 



s IJ = W^ J ] = { [LU) « 

where 

(L U )a B = i (■ 



B 

A 





A 

B' 



^JA'B _ n J ^IA'B 
AA' a AA'° 



(R") A B' = i (° IA ' A ° J AB> ~ ° JA ' A ° l AB, 



In this way the Dirac 4x4 representation decomposes into a left- and right-handed representation. Since 
primed and unprimed indices arc different kinds of indices the ordering does not matter. However, if we want 
the spinors a 1 AA , and g ja A to be the usual Pauli matrices it is necessary to have the primed/unprimed 
index as a row/column for g ja and vice versa for a 1 AA , [ill. Furthermore, g 1 AA , and <r are in fact the 
same spinor object if we use e A B and tA'B' to raise and lower the indices. Nevertheless, it is convenient for 
our purposes to keep the bar since that allows for a compact index-free notation g 1 = (1, g 1 ), a 1 — (1, — g 1 ), 
where the a 1 are (in matrix form) the usual Pauli matrices. 

The Dirac spinor can now be understood as a composite object: 

* = (£') (A.1) 

where §a and x A are left- an d right-handed spinors respectively. In this paper we take the left-handed 
component as encoding the quantum state. However, we could equally well have worked with the right- 
handed component as the result turns out to be the same. 

Although cab and Ea'B' are the only invariant objects under the actions of the group SX(2,C), the hybrid 
object a IA A plays a distinguished role because it is invariant under the combined actions of the spin-1 and 
spin-i Lorentz transformations, that is 

-I A' A , a J a A T A' ^JB'B ^IA'A 

a —f A j A B A qiG — G 

where A 1 j is an arbitrary Lorentz transformation and and A A B , are the corresponding spin-| Lorentz 

boosts. A^ is the left-handed and A B A (— KT 1A B ,) the right-handed representation of SX(2,C). 

The connection between 50(1, 3) vectors in spacetime and SL(2, C) spinors is established with the linear 
map g iaa , a hybrid object with both spinor and tetrad indices [53]. The relation between a spacetime 
vector (j) 1 and a spinor <f>A is given by 

= G IA ' A 4> A ,cb A . 

This relation can be thought of as the spacetime extension of the relation between 50(3) vectors and SU (2) 
spinors, i.e. this object is the Bloch 4- vector. This is in fact a null vector, and we can say that g ia a provides 
a map from the spinor space to the future null light cone. 
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A. 3 The geometric structure of the inner product 

In order to turn the complex vector space W into a proper Hilbert space we need to introduce a positive 
definite sesquilinear inner product. A sesquilinear form is linear in the second argument, antilinear in the 
first, and takes two complex vectors 4>a^a G W as arguments. The antilinearity in the first argument 
means that 4> A must come with a complex conjugation and the linearity in the second argument means that 
ipA comes without complex conjugation. In order to produce a complex number we now have to sum over 
the indices. So we should have something looking like (<fi\ip) = ^2 ^A'lpA- However, we are not allowed to 
carry out this summation: both the indices appear as subscripts and in addition one comes primed and the 
other unprimed. The only way to get around this is to introduce some geometric object with index structure 
jA'A e j^r* ^ yy* _ The inner product then becomes 

= I A ' A Ma- 

In order to guarantee positive definiteness, the inner product structure I A should have only positive 
eigenvalues. 

Now that we have defined an inner product structure we can state how the spinor index notation is 
related to Dirac bra-ket notation used in standard quantum theory. We can readily make the identifications 

|0) -0.4 (t\~i A ' A 4> A ,. 

In non-relativistic quantum theory, one would choose the inner product as I A A — S A A where S A A is the 
Kronecker delta. However, a different structure arises from the inner product of the Dirac field in the WKB 
limit. To see this we begin with the conserved current j M = ^'(x)j ,i '^>(x). The net 'flow' of this current 
through an arbitrary hypersurface forms the Dirac inner product, and is a conserved quantity. Now consider 



the Dirac inner product between two 4-spinor fields ^i(x), ^(aO in the Weyl representation ( A.l I. We have 
^>i(x)r^2(x) d£ M = / d» A , A x£(x)x£\x) + ^ AA '^ B ,(x)<l> 2 A (x) dS M (A.2) 



where the integration is over an arbitrary spacclikc hypersurface E. If n M is the unit vector field normal to 



the hypersurface and dS is the induced volume element, we write dS M = n M dS. (A.2) is further simplified 
by making use of the equations of motion mx A — ia^ A A D^(j)A, where the covariant derivative reduces to 
D^A ~ kfj,(f>A in the WKB approximation. In this approximation we obtain 

^(x)r^2(x) dS p a J uiula aB ' A ^ AA ,^ A ' B ^,{x)<t>l{x) + a» AA ' ^ Al (x)<f>l(x) dS M . (A.3) 

If we further assume that k\ = k 2 al i.e. the 4-momentum of the fields and ^2 in the WKB limit coincide, 
the inner product can be further simplified to 



9i(xypMx) dS M = J HZ ^OrVKzK d£ M (A.4) 
where we have made use of the identity Eqn (2.52) pl6] 

-aB'A h -BA'B au-BB'B aB-uB'B . Bu-aB'B . • aaBj - B' B f\ r\ 

(7 a AA' a — 9 a ~ 9 a + 9 a + 16 <7 7 • (A.OJ 

We therefore see that the inner product for the Weyl 2-spinor I A A , which we obtained in the WKB approx- 
imation, naturally emerges from the inner product of the Dirac field as the object contracting the spinor 
indices A', A at each point x of the fields. W\ 



B Jerk and non-geodesic motion 

We have seen that the transport of qubits as massive fermions is governed by the spin-^ Fermi- Walker 
transport equation (?J5|. One might then expect that the transport of qubits as polarization of photons along 

23 Note that the factor of 2 arises from differences in defining normalization: the Dirac spinor is normalized by ^(x)^(x) = 1 
with = (0, x), but the Weyl 2-spinor is normalized by o(x) &(x) = 1. 
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non-geodesic null-trajectories should similarly be governed by a kind of Fermi- Walker transport. Transport 
of polarization vectors for these non-geodesic null trajectories was developed by Castagnino |86[ 1871 188) . 
However, these two proposals are mathematically distinct, and it is not clear to us which one is the correct 
one. Furthermore, both of these proposals involve the 'jerk' along the path, i.e. the time derivative of the 
acceleration, making the transport equation for non-geodesic paths looks rather unpleasant. 

It is easy to show that any transport of a polarization vector along a null path must involve three or 
more derivatives of the trajectory ^(A), i.e. involve one or more derivative of the acceleration a M (A). From 
linearity and the requirement that the transport reduces to the parallel transport for geodesies we deduce 
that the transport must have the form — DX V — + T j<f> = 0. We now show that no such choice of T j 
containing only the 4-velocity and the acceleration exists that preserves the orthogonality <j> uj = between 
the 4-velocity and the polarization vector. We have that 

= cT\^ /w/) = = ^\ UI + ^ ai = ~ tI j Ui ^ J + ^ ai - 

However, if we now assume that the transport contains at most the second derivative of x^(A) (i.e. the 
velocity u 1 and the acceleration a 1 ) we deduce that T I J — a,u l uj + fiu 1 a,j + ^afuj + 5a 1 a,j. But since 
u 1 ai e we see that T I J ui<j) J = and we have thus deduced that (j) 1 aj — for all trajectories and all 
polarization vectors which is false. Therefore, we have a contradiction and we have to conclude that T I J 
contains one or more derivatives of the acceleration a 1 . 

It is, however, not clear that it is appropriate to study transport of polarization vectors along non-geodesic 
null trajectories. Physically, non-geodesic paths of photons can only be achieved in the presence of a medium, 
in which case the photon trajectories will be timelike. In our approach a physically motivated way to obtain 
non-geodesic trajectories would be to introduce a medium in Maxwell's equations through which the photon 
propagates. Nevertheless, even without explicitly including a medium, it is easy to include optical elements 
such as mirrors, prisms, and other unitary transformations as long as their effect on polarization can be 
considered separately to the effect of transport through curved spacetime. 
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